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Abstract. We review the density of states and related quantities of quasi one-dimensional 
disordered Peierls systems in which fluctuation effects of a backscattering potential play 
a crucial role. The low-energy behavior of non-interacting fermions which are subject to 
a static random backscattering potential will be described by the fluctuating gap model 
(FGM). Recently, the FGM has also been used to explain the pseudogap phenomenon in 
high-Tc superconductors. After an elementary introduction to the FGM in the context of 
commensurate and incommensurate Peierls chains, we develop a non-perturbative method 
which allows for a simultaneous calculation of the density of states (DOS) and the inverse 
localization length. First, we recover all known results in the limits of zero and infinite 
correlation lengths of the random potential. Then, we attack the problem of finite correlation 
lengths. While a complex order parameter, which describes incommensurate Peierls chains, 
leads to a suppression of the DOS, i.e. a pseudogap, the DOS exhibits a singularity at the 
Fermi energy if the order parameter is real and therefore refers to a commensurate system. 
We confirm these results by calculating the DOS and the inverse localization length for finite 
correlation lengths and Gaussian statistics of the backscattering potential with unprecedented 
accuracy numerically. Finally, we consider the case of classical phase fluctuations which apply 
to low temperatures where amplitude fluctuations are frozen out. In this physically important 
regime, which is also characterized by flnite correlation lengths, we present analytic results 
for the DOS, the inverse localization length, the speciflc heat, and the Pauli susceptibility. 
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1 Introduction 

As the temperature is lowered, some inorganic and organic conductors with a highly 
anisotropic crystal and electronic structure become unstable and undergo a Peierls 
transition, i.e. they develop a charge-density wave. This instability is due to their 
quasi one-dimensional nature which results in a (perfectly) nested Fermi surface. A 
qualitative understanding of the Peierls instability can already be gained by treating 
the phonon field of a quasi one-dimensional electron-phonon system in a mean-field 
picture ||l|-p. However, because of reduced dimensionality, fluetuations of the phonon 
field which can be identified as the order parameter field A(x) are crucial and signifi- 
cant deviations are to be expected. 

In a seminal paper, Lee, Rice and Anderson |9| introduced the one-dimensional 
so-called fiuctuating gap model (FGM), in which fluctuations of the phonon field are 
modeled by a static disorder potential. Calculating the leading-order correction of the 
electronic self energy of an incommensurate chain which is described by a complex 
order parameter field with (A(.t)) = and {A{x)A* {x')) — e'^^~^ 1/^, where f is 
the temperature-dependent correlation length, Lee, Rice, and Anderson obtained an 
approximate expression for the density of states (DOS), showing a suppression of the 
DOS near the Fermi energy, which is called a pseudogap. 

A few years later, Sadovskii apparently obtained an exact expression for the 
Green function of the FGM using Gaussian statistics for the higher correlation func- 
tions of the order parameter field which he could assume to be real or complex, referring 
to a band filling being commensurate or incommensurate with the underlying lattice. 
Recently, the experimental observation of a pseudo-gap state in the overdoped euprates 
above the superconducting phase transition led to a reincarnation of the FGM and 
Sadovskii's exact solution in the field of high-temperature superconductivity p]-|l3|. 
However, the revived interest in Sadovskii's solution also brought to light a subtle 
error in this solution which questions not only the solution itself, but also the 
work based on it. 
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Besides the limit ^ ^ oo where Sadovskii's solution is indeed exact ||lJ,|T^], Sadov- 
skii's solution can also be easily tested in the white-noise limit ^ — > 0, keeping D = A^^ 
constant, such that (A(a;)A* (x')) = 2DS{x — x'). Solving a stationary Fokker-Planck 
equation, Ovchinnikov and Erikhman 16 1 obtained an exact expression for the DOS 
for real A(a::). They showed that for small lo and (A(a;)) = 0, the DOS diverges 
as {p{yj)) oc |a;ln'' Singularities of this type at the band center of a random 

Hamiltonian have been discovered by Dyson in the fifties and have recently also 
been found in one-dimensional spin-gap systems 19|. It is important to note that 
in the FGM, the singularity is a consequence of phase resonance, and is not related 
to concrete probability properties of A(x) [^|2l|. In particular, the singularity is not 
an artifact of the exactly solvable limit ^ considered in Ref. ||l^. As argued in 
Ref. it is therefore reasonable to expect that for any ^ < oo the average DOS of 
the FGM exhibits a singularity at a; = 0. This general argument is in disagreement 
with Sadovskii's solution which for large but finite ^ shows a pseudogap and no 
singularity. In this work we shall reexamine the DOS of the FGM which determines 
the whole thermodynamics of the FGM and resolve the above contradictions. This 
emerged from the PhD thesis of the author |Q . Some of the results presented here 
have been published in a series of recent research articles [p2lp^,[2^. 

The organisation of this article is as follows: In Section we give an elementary 
introduction to diordered Peierls systems whose Fermi wave vector can be commensu- 
rate or incommensurate with the underlying lattice structure. Starting from a Frohlich 
Hamiltonian which describes a one-dimensional electron-phonon system, we will in- 
troduce the fluctuating gap model (FGM) as a low energy model in which the phonon 
system is essentially replaced by a fluctuating static backscattering potential which 
will serve as the order parameter field. 

Section || focuses on the one-particle Green function of the FGM. After calculating 
the Green function in the leading-order Born approximation which reproduces the 
result obtained by Lee, Rice and Anderson, we will develop a formally exact non- 
perturbative expression of the Green function as a functional of the disorder potentials 
based on a non-Abelian generalization of the Schwinger-ansatz. To calculate the DOS 
and inverse localization length, the introduction of phase variables will turn out to be 
very convenient. While one phase variable is simply related to the integrated DOS and 
satisfies a non-linear equation of motion which is equivalent to a Riccati equation, the 
other phase variable is related to the inverse localization length and can be expressed 
in terms of the first phase variable. These equations of motions will serve as the 
starting point for detailed calculations of the DOS and inverse localization length for 
various probability distributions of the disorder potentials in the next sections. 

In Section U, we will review known exact results of the DOS and inverse local- 
ization length in the limit of infinite correlation lengths and in the white noise limit. 
Generalizing the phase formalism developed in Ref. ||2^ such that A(a::) is allowed to 
be complex, we will derive a linear fourth-order Fokker-Planck equation previously 
only obtained within the framework of the method of supersymmetry |p6|| . The solu- 
tion of this stationary Fokker-Planck equation encapsulates all known results for the 
DOS and inverse localization length of the FGM in the white noise limit including the 
above mentioned Dyson singularity in the Ovchinnikov and Erikhman limit. Results 
for the case of infinite correlation lengths will finally be obtained by averaging the 
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DOS and the inverse localization length calculated for a constant disorder potential 
over an appropriate probability distribution of the disorder potentials. 

The case of finite correlation lengths of the order parameter field will be attacked 
in Section |^. Considering the equation of motion related to the integrated DOS, we 
will first argue that we expect for any finite ^ a Dyson singularity in the DOS. We 
will then set up an algorithm based on the equations of motion derived in Section 
^ which will allow for a simultaneous numerical calculation of the DOS and inverse 
localization length for arbitrary disorder potentials with unprecedented accuracy. For 
complex A.{x), Sadovskii's solution is not too far off from our numerical solution. 
In particular, for large correlation lengths A^^ ^ 1, the DOS at the Fermi energy 
vanishes as p{0) oc (A^^)"*'-^'* instead of p{0) oc (A^^)"^/^, as predicted by Sadovskii. 
However, for real A, we will find a pseudogap in the DOS for A^^ 3> 1 which for any 
finite ^ is overshadowed by a Dyson singularity of the form p{uj) ~ A jwln"^^ l"^!!^^! 
where A and the exponent a depend on the correlation length ^. As the correlation 
length ^ increases, a assumes the finite value a = 0.41, but the weight of the Dyson 
singularity vanishes with increasing correlation length. At the end of Section ^, we 
shall also discuss the case of only phase fluctuations of the order parameter which 
applies to sufficiently low temperatures where the amplitude of the order parameter is 
confined to a narrow region around A^ such that Gaussian statistics do not apply any 
more. We will find exact analytic results for the DOS and inverse localization length 
and we will also calculate the low-temperature Pauli paramagnetic susceptibility and 
the electronic low-temperature specific heat. 



2 Peierls systems and the fluctuating gap model 

In this introductory section we will replace the dynamic phonon ensemble in Peierls 
systems by a static backscattering potential which we will identitfy as the order param- 
eter. Due to reduced dimensionality, fluctuations of this order parameter field are very 
important and we will determine its statistics in the context of a brief discussion of 
a generalized Ginzburg- Landau functional. Finally, we introduce the fluctuating gap 
model (FGM) as a low-energy model which takes into account these fluctuations. 



2. 1 Frdhlich Hamiltonian and Peierls instability 

The formation of periodic lattice distortions and charge-density waves in Peierls chains 
is due to the electron-phonon interaction in these quasi one-dimensional materials |l|- 
ill^Hll- Since particle- hole excitations with momenturrj^ 2kF are possible for very 
small excitation energies, the Lindhard density-density response function exhibits a 
singularity at q = 2kp. Kohn showed that this singularity should be conveyed into 
a kink in the phonon spectrum pTj . While these Kohn anomalies are rather weak in 
isotropic materials, they can lead to a substantial alteration to the phonon dispersion 
in quasi one-dimensional materials with a topology of the Fermi surface which shows 
perfect nesting. At low enough temperatures, the renormalized phonon mode at 2kp 
can scale all the way down to zero, i.e. become gapless. This process is called softening 



In this work wc choose units such that h - 
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of the phonon mode. Since Wrcn(2fcF) ^ 0, a static lattice distortion with wave vector 
2kF may now arise. Simultaneously, there is a formation of a charge density wave. As a 
consequence, the discrete translational invariance is broken. The same physics can also 
be described by considering the thermodynamics of a Peierls system. This approach 
will also allow to go beyond a mean-field picture and will therefore be followed here. 

A Hamiltonian to describe a one-dimensional electron-phonon system was proposed 
in 1954 by Frohlich |§: 

^ = E + E ^« + T.^s\{h, + bl^) . (1) 

k,(j q g 

The system has length L = Na where a is the lattice spacing, and periodic boundary 
conditions are assumed. c\, and Ck are fermionic creation and annihilation operators 
with momentum fc, spin cr, and energy efe. While Ck = /2m for free electrons, in the 
tight-binding approximation one has Cfc = —2tcoska. The second term in the Frohlich 
Hamiltonian (]l) describes phonons with phonon dispersion ujq. 5j and bq are bosonic 
creation and annihilation operators with momentum q which is confined to the first 
Brillouin zone. For a chain with only nearest-neighbor interactions we have (see for 
example Ref. |^^) cuq = 2a;o sin |(jra/2|. Finally, the last term in Eq. (|^) models the 
interaction of the phonon system with the fermions. The phonons are linearly coupled 
via the electron-phonon coupling constant gq to the Fourier components of the electron 
density 

^9 = E 4+9,<TCfc,^ • (2) 

The phonon operators bq and 6j are directly related to the operators of the normal 
coordinates Uq of the lattice system by 



1 



1/2 



bg + blq] . (3) 



2MLUqJ V 

Here, M is the ionic mass. The lattice displacement operators of the ions at x„ = na 
are given by its Fourier transform, 

^M=Y.^^^-'{^Y\bq+blq) . (4) 



9 



In a mean-field picture, the Peierls transition will lead to a non-vanishing expectation 
value (u(x„)) which implies that the system exhibits a static lattice distortion P, p3[. 
This static lattice distortion is accompanied by a single-particle gap, a charge density 
wave and unusual electronic transport properties ^ . 



2.2 Euclidean action 

In an Euclidean functional integral approach, the Frohlich Hamiltonian is conveyed 
into the action (see, for example, Negele and Orland [pl[) 
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where 




(8) 



Here, /3 = l/T is the inverse temperature and efc = — /i is the energy dispersion 
measured with respect to the chemical potential /i. While the conjugated Grassmann 
variables '0fc,ii'„ V'fc li describe fermions with momentum k and fermionic Matsub- 
ara frequency (D„ = (2n+ l)7r//3, bq^^j^ and 6* ^^^ are complex (bosonic) phonon fields 
with momentum q and bosonic Matsubara frequency ujm = 27rm//3. In terms of the 
above action, the partition function reads 

V{r,i^}'D{b*,b}exp[-S{r,i';b*,b}] , (9) 

where V {tp* , ^} and V {b* , 6} are appropriately normalized fermionic and bosonic inte- 
gration measures ||3l) . Using the variable transformation = {pq,i^m +^-q.-w„Ji 
riq.u,^ = {bq.^u^m such that (j)*^^^^ = ^-g.-c^™ and r/*^„^^ V-g-uj^, 
7] may easily be integrated out resulting in 

5{0*,^;0} = 5ci{V*,V'} + VW + ^int{V'*,V';0} , (10) 
where Sc\{ip^^ip} is unchanged and 

1 1 (jj^ ~\~ up" 

V'; 0} = E I E V':+,,^„+.„V'fe,i.„ I ■ (12) 

So far, no approximation has been made. In the following, we will restrict ourselves to 
the low-energy physics of the weak-coupling limit, so that only fermions in the vicinity 
of the Fermi energy are involved. In this case the Fermi energy may be linearized 
around the two Fermi points, such that it assumes the form 

lk = VF{\k\-kp) . (13) 
To separate right- and left-moving Fermions, let us introduce the spinor field 

^ ^ { \ ^ { \ (14) 
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and its conjugated counterpart 

The electronic part of the action may easily be rewritten in terms of these spinor fields 
and the inverse non-interacting Matsubara Green function 

iun — Vpk 

iuJn + Vpk 

Since the momentum transfer of the phonons is either small compared with the Fermi 
momentum or approximately 2fci?, we decompose 4>q.uj„^ according to 



Go ^(^''^«) = 



(16) 



V 



Va. 



directly translates into V*^^ 



t>q-2kF,iAj„ 



(17) 



such that \q\ < kp. While 4>q^^^ = 4>-q,-uj„ 

a similar relation for A^^^^^ does only hold ii'ikp is a reciprocal lattice vector. A* — 
A_q is therefore only true for a half- filled band for which -nja — 2kp. We will 
refer to this case as the commensurate case. The more general case for which kpa/n 
is an other fractional number is also called commensurate but will not be discussed 
here. In the incommensurate case, for which kpa/ir is well separated from any simple 
fractional number, all A 



and A_g',_c 



are independent. We will see in this 
work that commensurate and incommensurate Peierls systems can have very different 
physical properties. 

Defining the matrices Gq ^ and V via 



G 







'k.k' 



' k.k' ,LUn -^n' 

our action turns into 



,G„ ^{k,uj„ 



S{^\t, V, A,A*} = 5ei-phW\ V'; V, A,A*} + ^ph{V, A, A*} 



where 



k,k' ,UJn :^n' '"^ 

SMV.A,A*}^lpLj:J-^ 

1 



k.k' 



Ipk 



(18) 
(19) 

(20) 

(21) 



q,UJ-n 

C ^-^ 



V* V 



\g2kF+q\ 



'^2kF+q 



^2kF+q 



A* 

9,^77 



A 



and 



commensurate case (half-filled band) , 
incommensurate case . 



(22) 



(23) 



While in the incommensurate case ±2kp lie (up to a reciprocal lattice vector) inside 
the first Brillouin zone, ±2kp lie directly on the border of the first Brillouin zone in 
the commensurate case. In this case the factor of 1/2 in the last line in Eq. ( p2| ) avoids 
overcounting. 
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2.3 Generalized Ginzburg-Landau functional 



Since the action S{-4>^ j^piV, A, A*} is only Gaussian in both the Fermion and the 
phonon fields, either of them can easily be integrated out. To derive a time-independent 
(generalized) Ginzburg-Landau theory which will enable us to determine the statistics 
of the phonon field, we ignore quantum fluctuations, i.e. we ignore all terms involving 
finite bosonic frequencies. Because we will be only interested in static properties of the 
Peierls system, this should be a reasonable approximation for not too small tempera- 
tures. Integrating out the fcrmionic degrees of freedom, the action -0; V, A, A*} 
turns into the free energy functional f3F{V, A, A*}. 

For the derivation of a generalized Ginzburg-Landau functional in which the free 
energy functional is expressed in terms of gradients of the order parameter field 



A(x)=E^'''^?' 



(24) 



we refer the reader to the PhD thesis of the author |2^ or to Refs. |32|-p8[ where the 
gradient-expansion is developed in the context of superconductivity. Up to terms of 
second order in the gradient, we find 



F{A, A*} = F("){A, a*} + F(2){A, a*} 



f(°){a,a*} 



spo 



dx 



F(2){A,A*} = 5po 



dx ■ 



27r 
T 



27r 
T 

E 

ii„>0 



0<a;„<eo 

1 |9xA| 



'^'^.-t-iAr 



2A 
lAPl- 



{^i + \A\^r 



32 



{^i + \A\^r 



(25) 
(26) 

(27) 



For i^(°){A, A*} to be finite and to avoid logarithmic divergences, the sum in Eq. ( [2^ ) 
needs to be regularized by an ultraviolet cutoff eg. Although the ultraviolet cutoff eo 
was introduced here in the sum over Matsubara frequencies instead of as a cutoff in 
the momentum integral, expanding Eqs. (^6|) and (|2^) in the regime /3|A| ^ 1, we get 
the usual Ginzburg-Landau functional, 



^{A,A*}.ffI 



dx 



a{T) \A{x)\' + b{T) |A(x)r + c(T) \d,A{x)Y 



where the coefficients a(T), h{T) and c{T) are given by 
T 



a{T) 

b{T) 

c{T) 
Here, C(3) ; 



In 



rpMF ' 



T, 



MF 



1.134eoexp(-l/A) , 



/ 1 



V47rT 



7C(3) , 



A=- 



\A-kTI 

1.2 is the Riemann zeta- function of 3 and 

2 



(28) 

(29) 
(30) 
(31) 

(32) 
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is the dimensionless coupling constant. is the critical mean- field temperature for 
the Peierls distortion. If we include spin, s = 2, otherwise we have s = 1. 

2.4 Breakdown of the mean-Held picture 

The experimentally observed Peierls transition, including the static lattice distortion, 
the charge-density wave and the occurrence of a single-particle gap can qualitatively 
already be understood in a mean-field picture. However, due to reduced dimension- 
ality, fluctuations of the order parameter field are very important and dramatically 
change this scenario. The Mermin- Wagner theorem states that these fluctuations lead 
to the absence of long-range order, even at very low temperatures This pre- 

cludes a spontaneously broken continuous symmetry. But how can one explain the 
experimentally observed charge-density wave which breaks a continuous translational 
symmetry in strongly anisotropic materials like blue bronze ^? The answer is simply 
this: These materials are quasi one-dimensional, but not strictly one-dimensional. As 
we will see below, in a strictly one-dimensional material, the correlation length S,{T) in- 
creases with decreasing temperature, but for any finite temperature cannot approach 
infinity. At very low temperatures, however, even a very weak interchain-coupling 
can lead to the onset of three-dimensional order such that the system can undergo a 
Peierls transition. Of course, the transition temperature is not the mean-field transi- 
tion temperature . Lee, Rice and Anderson pointed out that one should expect 
T^^ « jT^^ . For a derivation of an adequate three-dimensional microscopic theory 
see McKenzie |Q and references given therein. Here, we will especially be interested 
in the temperature regime above the Peierls transition so that it suffices to consider 
only the strictly one-dimensional case. 

2.4.1 Correlation functions of the order parameter field 

We will now consider the fluctuations of the order parameter field A(x) and calculate 
the correlation functions of A(a::) which describe the phonon statistics. 

Gaussian approximation 

At temperatures far above the mean-field critical temperature T^^, the coefficient 
a(T) becomes large enough such that anharmonic corrections to the free energy func- 
tional may be neglected. Truncating the free energy functional ( p8| ) at the second 
order, we find 



where A^^{T) = and (-^{T) = (^) . While (A(x)A(x')) vanishes 



for complex A, it is equal to {A{x)A*{x')) for real A. In the considered Gaussian 
approximation, higher correlation functions are simply given by Wick's theorem. The 
Gaussian approximation is good at sufhciently high temperatures and is the usual 
approximation made when higher correlation functions are too complicated or not 
known. 



{A{x)) = , 
{A{x)A*{x')) = A2(r)e-l^-^'l/«(^) 



(33) 
(34) 
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Anharmonic corrections and the case of only phase fluctuations 

As the temperature is lowered and approaches the mean-field critical temperature 
T^^, fluctuation effects become important and the mean-field picture breaks down. 
However, as shown by Scalapino, Sears and Ferrel p| using the transfer matrix tech- 
nique, the first two moments of a one-dimensional Ginzburg-Landau theory are still 
approximately given by Eqs. ( |33| ) and (|3^). For temperatures well below T^^, one 
finds for real A{x) an exponential increase of the correlation length with decreasing 
temperature, while for complex A(x) the correlation length increases as the inverse 
temperature. This last result can be understood as follows: For small temperatures 
T <C , the generalized Ginzburg-Landau functional has the shape of a "Mexican 
hat" and is dominated by its minima. Amplitude fluctuations get gradually frozen out 
and, for complex A(x), the phase of the order parameter A(x) w A^e*''*^^^ fluctuates 
at the bottom of the "Mexican hat" . Ignoring, as before, quartic terms in the gradient 
expansion of the free energy, the free energy is given up to an irrelevant constant by 

^(phasc)|^| ^ FiP^^^'){d,§/2} = I spsiT) [ dxV^ix) . (35) 

2 Jo 

In analogy to the theory of superconductivity, 

V{x) = d^d{x)/2 (36) 
can be interpreted (up to a constant 1/m*) as the superfluid velocity and 

Ps{T) = p,'^-^ ^ (37) 

P ^t^oR + IA.P)^ 

can be interpreted as the superfluid density. Formally, the free energy is identical to 
the kinetic energy of a superflow. A two-dimensional analogue of Eq. ( |35| ) has been 
used by Emery and Kivelson 43 1 in their theory describing superconductors with 



a small phase-stiffness. At T = 0, the sum in Eq. (37) turns into an integral which 
can be done analytically and gives /9s(0) = po, i-c. at T = 0, the superfluid density is 
equal to the density of states. Plots of Ps{T) for As{T) given by the BCS gap equation 

1 o„ 1 

(38) 



1 _ 27r ^ 1 



and for A^ = 1.76r^^^ are shown in Fig. |l|. 

Since F^^^'^^'^^ {V} is only quadratic in V{x), correlation functions of A(a;) can 
easily be calculated. The first two moments are given by Eqs. (|3^) and (Q), where 

(39) 

For T < T^^/A we have psiT) sa p^iO) = I/tt, such that in this strictly one- 
dimensional theory we find ^(T) — s/2tiT cx 1/T which for fermions with spin 1/2 
agrees with Griiner's lil result ^(T) = I/ttT. 
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Fig. 1 Plot of the superfluid 
density as a function of temper- 
ature for As = 1.76T^^ (dashed 
line) and As{T) given by the 
BCS gap equation (^) (solid 
line) . 



2.5 The Hamiltonian of the Huctuating gap model 

The fluctuating gap model (FGM) describes electrons subject to a static disorder po- 
tential which can be seen as an approximation of a phonon field. For a particular 
realization of the disorder, the electronic part of the action given in Eq. (|2|) corre- 
sponds to the Hamiltonian 



fe.fe' 



where 



Hi. 



k,k' 



vpkSk.k' + Vk-k> 
A* 



C+.k' 



-(k~k') 

A Fourier transformation leads to 



^k-k' 

-vpkShk' + Vk-k' 



n 







Tp+ix) 
V'- {x) 



with 



H{x, —idx) — —ivpOxCs + V{x)aQ + A{x)a^ + A* (a;)cr- 



(40) 



(41) 



(42) 



(43) 



This is the Hamiltonian of the FGM. cTi are the usual Pauli matrices, (Tq is the 2x2 unit 
matrix, and a± = ^(tTi ± ia2). Recall that we have linearized the energy dispersion 
such that the FGM can only describe the low-energy physics of Peierls chains in the 
weak-coupling regime. As a further approximation, we have considered the phonon 
field to be static. It will now be our aim to calculate disorder-averaged quantities for 
the model described by this Hamiltonian. As we will discuss in Section ^, instead of 
averaging over the disorder, it is also possible to consider a typical realization of the 
disorder potential. 

2.5.1 The fluctuating gap model in other physical contexts 

In this section, the fluctuating gap model (FGM) emerged as an effective low-energy 
model to describe quasi one-dimensional materials which undergo a Peierls transition. 
Our strictly one-dimensional theory applies to temperatures above the Peierls tran- 
sition before three-dimensional fluctuations become important and eventually lead to 
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a phase transition. Formally, the Hamiltonian of the FGM is of the Dirac type and 
describes electrons in a disordered potential. In the language of relativistic quan- 
tum field theory, the backscattering potential can be interpreted as a random (and 
complex) mass. The FGM has also applications in other fields of physics. As shown 
in Refs. the Hamiltonian of disordered spin Peierls systems [p|, p8|-|50| can 

be mapped by a Jordan- Wigner transformation onto the Hamiltonian of the FGM. 
In a semiclassical approximation of superconductivity, it is also possible to replace 
the original three-dimensional problem by a directional average over effectively one- 
dimensional problems |^ which in the weak coupling limit are described by the FGM. 
This method has been used in Refs. |36-38| to derive the gradient expansion of a clean 
superconductor. A generalization of the FGM towards higher dimensions to describe 
the phase above the phase-transition in underdoped high- Tc superconductors by anti- 
ferromagnetic short-range order fluctuations was considered in Refs. |ll|-|l3|. 



3 The Green function of the fluctuating gap model and related quantities 

In this section, we will introduce different concepts to calculate the Green function and 
related quantities of the fluctuating gap model. The density of states and the localiza- 
tion length will be of special interest. In particular, we will develop a non-perturbative 
method which allows to calculate these quantities simultaneously for arbitrary given 
disorder potentials. 

3.1 The retarded Green function 

In the following, we are going to consider the retarded Green function Q^{x,x';u!) of 
the fluctuating gap model. This retarded Green function is of special interest because 
it can be related to several quantities which are in principle experimentally accessible. 
The trace of the imaginary part of the Green function at coinciding space points 
determines the local density of states (DOS), 

p{x, Lu) = -7r-ilmTr[^;^(a;, x; lo)] . (44) 

Averaging p{x,Lu) over all space points gives the DOS p{lo), which is the fundamental 
quantity that determines the whole thermodynamics of the FGM. It will turn out that 
the trace of the energy-integrated space averaged Green function at coinciding space 
points T(uj) will be easier to calculate than its non-integrated form. While its imag- 
inary part is proportional to the integrated DOS A/'(w), the Thouless formula states 
that Rer(cj) is equal to the inverse localization length £~^{uj). As in the following 
sections, we will usually only consider the DOS and the inverse localization length for 
positive frequencies w. Due to particle-hole symmetry, the DOS and the localization 
length are symmetric with respect to the Fermi energy so that after setting the Fermi 
energy equal to zero we have p{u!) = p{—oj) and £~^{uj) = £^^{—u!). It therefore 
suffices to consider the case a; > 0. 

The retarded 2x2 matrix Green function Q^{x, x';uj) to the Schrodinger operator 
uj — H satisfies the differential equation 

[uj + iQ+ - H{x, -id^)] g^ix, x';io)= 5{x - x')aQ . (45) 
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The positive but infinitesimal imaginary part added to the frequency lo indicates that 
we have to impose the correct boundary conditions applying to a retarded Green 
function. 

3.1.1 Free fermions 

It is easy to calculate the Green function for free fermions. In this case, V{x) = 
A{x) = 0, such that the system is translational invariant, and Eq. (45) simplifies to0 

[lo + iO^ + ia-zdx] Goix — x';lu) — ao6{x — x') . (46) 

Taking the Fourier transform of this equation from real space to momentum space 
gives 

[w + iO+ - kas] Q^ik; uj) ^ . (47) 

Q^{k]Ljj) = J dx e~'^''^Go'{x;uj) can now be found by a simple matrix inversion. If 
a — I accounts for right- and a = — 1 for left-moving fermions, the matrix elements 
of GQ {k;uj) are given by 

(eo«)_,(fc;.) = (48) 

A simple Fourier transformation back to real space now gives the retarded propagator 
of free fermions in real space, 

z(c;o«)_,(x;c.)=<5„^„,0(ax)e-"-. (49) 

Here, d(x) is the Heaviside step function 

oi^)^{l ; m . (50) 

For concreteness, let us also define 9{0) = \imx^o[0{x) + d{—x)]/2 = 1/2. While the 
matrix elements at x = are sensitive to the definition 9{0) = 1/2 which amounts to 
defining ^^(a; = 0; w) = i lim^^o+ [G^'i+x; to) + G'^i-x; lo)] , the local DOS p(x, lo) = 
— 7r~^ImTr[C/(5^(0; [jj)] does not depend on this definition because it only involves the 
harmless quantity 9{x) +9{—x) = 1. Due to translational symmetry, the total DOS is 
equal to the space-independent local DOS, 

Po{lo) = TT-' . (51) 

Note that the DOS of free fermions is independent of the frequency because we have 
linearized the energy dispersion. 

A further Fourier transformation of Eq. ( |4^ ) from frequency to time gives 

* t) = So.,a'0mc^x - t) . (52) 

This free retarded Green function in space and time allows for a simple interpretation: 
A fermion put into the system at t' = as a right- or left-mover will at time t > have 
traveled a distance |x| = i = vpt in the positive or negative direction, respectively. 
The fermion can not be observed in the system at times t < 0. 

^Besides h and kg, we also set the Fermi velcity vp equal to one. 
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Fig. 2 Diagrammatic expansion of the matrix Green function. While the single line repre- 
sents the Green function of free fermions, the double line is a graphical representation of the 
full Green function Q^{x,x';lu). The crosses denote the disorder potential V(a;). 



3.2 Dyson equation and perturbation theory 

One way to handle the disorder is to consider the disorder potential as a perturbation 
and expand the Green function in powers of this potential. Defining V(a;) = V{x)aQ + 
A(a;)CT_|_ + A*(x)cr_, Eq. ( p5| ) may be written as 

[iasd^ +L0 + iO+] g"{x, x'] uj) = d{x - x') cto + V(a;) g'^ix, x'] to) . (53) 

Substituting x hy xi, multiplying the resulting equation from the left with the free 
Green function g^ix — xi\lo) and then integrating over xi gives the Dyson equation 

g^{x,x';uj)=g^{x-x';uj)+ [ dx g^{x- xi;cj)V{xi) g'^{xi,x';uj) . (54) 



Iterating this Dyson equation, the exact Green function can be expressed in terms of 
the free Green function and the disorder potential: 

oo 

g«(x,x';c.) = ^e^(x,x';a;), (55) 

n=0 

where g^^x^x' \lo) = g^ ix — x']uu) is the free Green function calculated above, and 
for n > 1 the fimctions gf^'{x,x';uj) are given by 

g^ix, x'; ^) = J dxi... J dxn 

g^(x - Xn, uj) V(x„) g^ixn - Xn-i;uj) ■ ■ ■ V(a;i) go{xi - x'\ uj) . (56) 

Recall that the right-hand side of this equation involves the product of 2 x 2-matrices. 
The perturbative expansion of the full Green function can be visualized by using 
Feynman diagrams (see Fig. ||). 

The physical interpretation of the perturbation expansion is simple: The pertur- 
bative expansion takes into account all possibilities of a particle moving through the 
sample getting scattered at the various static impurities. While A(a;) changes the 
direction in which the particle travels and therefore can be interpreted as a backscat- 
tering potential, V(x) does not change the direction of the particle, so that it only 
leads to forward scattering. 
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Fig. 3 Diagrammatic representation of the averaged matrix Green function. The single line 
represents the Green function of free fermions and in this case the (averaged) double line is 
a graphical representation of the full (averaged) Green function {Q'^{x,x' -jUj)) . The dashed 
line denotes the disorder average (V(3;)V(a;')). 

3.2.1 Boundary conditions of the retarded Green function 

Below, we will consider a non-perturbative approach to calculate the Green function of 
the FGM. The above perturbative expansion can be used to obtain the correct bound- 
ary conditions of the full retarded Green function: Let us consider {Gn )aa' (^'^'i'^)- 
According to Eq. (p6|), its expansion in a product of free Green functions and the 
static impurities starts with (Gq)^^ {x — Xn;uj) and ends with (Gq)^,^, {xi — x';lj). 
These terms are proportional to d{a{x — a;„)) and 9{a'{xi — x')), respectively, so that 
{^n)aa' (■^'■^'''^) ^'^ vanish as ax —oo or a'x' — > oo. Since this reasoning 
applies to all orders in perturbation theory, it also applies to the full Green function. 
If we demand the potentials to vanish outside the interval [—A, L + A], the boundary 
condition can also be written as 

{G^)^^,{-aA,x';iu)^0, (CJ^)^^, (x, a'(L + A); = . (57) 

3.2.2 Second order Born approximation 

Let us now consider the disorder-averaged Green function. As discussed in Section 
^, above the Peierls transition, the first two moments of the order parameter field 
A(a;) are given by (A(x)) = and (A(a;)A* (a;')) = A^e^l^'-^'l/^ . Following Lee, Rice 
and Anderson we ignore the forward scattering disorder, i.e. set V{x) = 0. For a 
perturbative approach we assume Gaussian statistics for the higher moments of the 
order parameter field such that these moments can be separated according to Wick's 
theorem. A diagrammatic representation of the averaged Green function is shown in 
Fig. 1^. An infinite number of diagrams can be summed up by introducing irreducible 
diagrams which by definition cannot be separated into two disconnected diagrams by 
cutting a single propagator. The corresponding amputated diagram is obtained by 
eliminating all outer propagators. The sum of all amputated irreducible diagrams is 
known as the self-energy and is diagrammatically presented in Fig. ^. In terms of the 
self-energy, the averaged Green function reads in momentum space 

(58) 



{G''{k;u;))= {G^ik;u;)) '^J:{k;u;) 



The simplest approximation to take into account fiuctuation effects of the order pa- 
rameter is to consider only the first diagram in Fig. H. This approximation is known 
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as the second order Born approximation and is essentially the approximation made 
by Lee, Rice and Anderson in their seminal paper |^ in which fluctuations of the 
order parameter of the FGM were taken into account for the first time. A special 
non-Gaussian probability distribution of A(a::) involving only phase fluctuations for 
which the second order Born approximation turns out to be exact is presented in Ref. 

Ill- 

3.2.3 The self-energy 

Since (V(a;)) — 0, the self-energy in the second order Born approximation is given by 
SsCa; - x'; uj) = {Y{x)g^ix - x'; c^)V(x')) . (59) 

Placing Eq. ( |49| ) into this equation, we get 

(Sb)„„, {x - x'- Lu) = e-l--'l/« {gR)_ _ ^) . 

= -iS^.o,, e{-a[x - a;'))e-»[-+V«](— ') _ (go) 

As one should expect, the process of averaging restored translational invariance. Tak- 
ing the Fourier transform of Eq. ( |60| ) , we arrive at 

(Sb)„„, {k-M = j dxe-^"^ (Sb),„, ix;u;) = 6^,^, ^^^^^^./^ ' (61) 

Within the second order Born approximation, we therefore flnd for the one-particle 
Green function 

■ (62) 

This result was first obtained by Lee, Rice and Anderson 

3.2.4 The density of states and the inverse localization length 
Integrating Eq. ( |6^ ) over k and taking the trace, we obtain 

TV g%{x, x; u) = -z^^±^^^^ , (63) 



+ 



+ 



+ • • • 



Fig. 4 Diagrammatic representation of the (irreducible) self-energy. As in the above Figure 
^, the single line represents the Green function of free fermions and the dashed line denotes 
the disorder average (V(a;)V(a;')). 
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Fig. 5 Plot of the DOS ps (i^^) and the inverse localization length Ig^ (lo) calculated for As^ — 
0.2, 0.5, 1.0, 2.0, 10.0, and oo (mean-field result) in the second order Born approximation. 



where -y/z is defined as the principal part of the square root with the cut chosen along 
the negative real axis. The imaginary part of Eq. (|6^) gives the (averaged) DOS, 



PB(tj) = pqRc 



V(c^ + z/2e)2-A2 



(64) 



As we will show in Subsection 3.5, the real part of the Green function is equal to 
the derivative of the inverse localization length l~^(u)). Integrating this equation with 
respect to and setting the integration constant at infinity equal to zero, we obtain 



1/2C 



(65) 



A plot of both pb{^) and 1'^^{lo) is shown for different values of the dimensionless 
parameter As^ in Fig. |^. With increasing correlation length ^, the DOS gets more 
and more suppressed for oj < A^. However, instead of a real gap the fluctuations can 
only create a pseudogap. At w = 0, the DOS is given by 



Pb(0) = po- 



1 



v/1 + (2A.C)^ 
such that for A^^ 3> 1, the DOS vanishes as 



Pb(0) 



1 



oc 



(66) 



(67) 



2A,^ A,$ 

As the correlation length approaches infinity, the DOS assumes the mean-field result 



A2 
^0 



At zero frequency, the inverse localization length is given by 
£b'(0) = A, V^1 + (1/2A,0'-1/2A4 



(68) 



(69) 
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Fig. 6 Plot of the spectral function psikp + k;uj) with Afc = 0.0 (left-hand side) and 
0.5 (right-hand side) calculated for As^ = 0.2,0.5,1.0,2.0,10.0 in the second order Born 
approximation. 



3.2.5 The spectral function 

Another interesting quantity related to the single-particle Green function is the spec- 
tral function 

piakF + k;cj) = -hm{g^ik;u;))^^^ . (70) 

Experimentally, the spectral function can be measured by angular resolved photoemis- 
sion spectroscopy (ARPES). It follows from Eq. (^) that in the Born approximation, 
the spectral function is given by 

PBiakp + k;uj) — pf) % . (71) 

{A1 - {uo^ - k^)f S,-^ + [uj - aky 

Plots of the spectral function pBiakp] and psiakp + k; lu) with k = 0.5As as func- 
tions of LU are shown for different values of Ag^ in Figs ^. While for small correlation 
lengths, i.e. A^i^ <^ 1 the spectral function exhibits a maximum near uj = ak, for large 
correlation lengths we find two maxima near uj — zt^/A^ + k^, the closest one to ak 
having the larger weight. 

Corrections to the second order Born approximation and Sadovskii's solution 
An attempt to sum up all diagrams in the perturbative expansion of the averaged 
Green function was made by Sadovskii in the late seventie s |10|] . The first diagram 



in Fig. ^ which represents the contribution given by Eq. (^) can be calculated by 
attributing a factor of A^ to the phonon line and a factor {Go')^^ {k;uj + i/C) to the 
fermion line. Sadovskii conjectured that all other non-vanishing diagrams could be 
calculated using a similar recipe. (In the general case, the imaginary part added to 
Lj has to be multiplied by the number of phonon lines above a given fermion line.) 
This recipe enabled Sadovskii to give a continued fraction representation of the self- 
energy (or Green function). Sadovskii's solution was known as the only available 
exact solution of the pseudo-gap state (see Ref. [Q) and was therefore also used 
by other authors ||ll], n2, bSl . However, only recently, Tchernyshyov |14| discovered 
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an unfortunate error in Sadovskii's solution which also turned the work based on it 
into question. As pointed out by Tchernyshyov, Sadovskii's conjecture already breaks 
down for terms contributing to the next to leading order in the self-energy. While 
Sadovskii neglects the second term in Fig. ^, which is correct for the incommensurate 
case, the last term in this figure already gives a different contribution than conjectured 
by Sadovskii. Instead of trying to correctly sum up all diagrams in the perturbative 
expansion of the Green function we will now develop a method which will allow for a 
non-perturbative calculation of the Green function. 



3.3 Non-Abelian Schwinger-ansatz 

We base our non-perturbative approach to calculate the Green function of the FGM on 
a matrix generalization of the Schwinger-ansatz . To make the differential operator 
—idx proportional to the unit matrix, we first factor out a Pauli matrix (73, so that 
the retarded Green function 

g'^ix,x';u;)^a3g'^{x,x';u;) (72) 

satisfies 

[idx " M{x, (jj + iQ+)]Q'^{x, x'; u) = S{x ~ x')ao , (73) 

where 

M{x,Lu) = [V{x) ^ uj + A{x)a+ + A* {x)a^]a3 (74) 
is a tracelcss matrix. We now try to solve Eq. (|7^) by making the ansatz 

g^{x,x',uj) = u{x,uj)g^ix~x')u-\x',uj) , (75) 

where U{x,uj) is an invertible 2x2 matrix and Gq {x) is the Green function to the 
operator idx + iO+(T3, i.e. 

The ansatz ( [T^ ) resembles the transformation law for the comparator in non-Abelian 
Gauge theory (see, for example, the book on quantum field theory by Peskin and 
Schroder |^ ) , and since it is also similar to the scalar Schwinger-ansatz ||5^ which is 
sometimes used in functional bosonization of interacting fermions , we will refer to 
it as the non-Abelian Schwinger-ansatz j22|,^,|3^. But note that in contrast to earlier 
formulations of the non-Abelian Schwinger-ansatz p2|,|3^, we use the zero- frequency 
free retarded Green function, such that the whole uj dependence is included in U{x,uj). 

In the following we are going to suppress the parameter w. The ansatz ( |75| ) indeed 
solves Eq. (|7^) if U (x) satisfies 

[idx - M{x)] U{x) = . (77) 

To establish our formalism, let us first restrict the disorder potentials to the interval 
(—A, L + A). While the potentials are assumed to be constant in the intervals (—A, 0) 
and (i, L + A), they are allowed to fluctuate in between. Later we can let A — > oo (or 
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we will set A = and let L — > oo). The boundary conditions for the retarded Green 
function given in Eq. (57) now renders into|^ 



?7i2(-A) = C/2i(i + A) = 



Two different solutions of Eq. ( 77 ) are given by 
U+{x) — Texp 

C/_(x) = T-^exp 



-I I M{y)dy 

/-A 



L+A 



M{y)dy 



(78) 

(79) 
(80) 



where Texp is the path-ordered and T ^exp is the anti-path-ordered exponential func- 
tion. Both, U-\-{x) and U-{x) can be expressed in terms of the ^'-matrix, 



S{x, x') = Texp 



M{y)dy 



(81) 



By definition, U+{x) = S(x,-K) and U-{x) = S^^{L + A, a;). Because = (J^May, 
and TrM = 0, the S'-matrices satisfy = cr^S^^a^ and dots' ~ 1, which means 
that they belong to the non-compact group SU{\, 1). It follows that the elements of 



S satisfy S22 = S^^^ S21 = 8^2, and \Sii 



\Si2\ = 1. While each Ua{x) only obeys 



one of the two conditions (iTq), the combination 



U{X) EE — 



U-nix) U+i2{x) 

U-2l{x) U+22{X) 



(82) 



satisfies both boundary conditions. Here, u — 6*22 (T+ A, — A) — J7_ii(— A) = U+22{L- 
A), so that detJ7(a;) — 1. Defining 



Uq = {Uall, ~Ua2l)'^ , Vq = (— C/q12 , C/a22) 

(such that = ctiU* ), we obtain from Eqs. ^72), 



5^) and m 



■rii/ ' \ ai ,^ ^-i^W+jx') , ^, , , -^+{x)xl{x') 
ly [x, X -jLj) ^ 0(x — X ) 1- v[x — X) 



(83) 



(84) 



and vl are the adjungated row vectors to the column vectors u+ and v_, so 
that u_u^ and v+vl are 2 x 2-matrices. Note that Eq. ( |8^ ) involves only U+12, 
U+22, UI12 and ?7^22j but not its complex conjugates. In principle, Eq. ( ^ ) allows to 
determine the full Green function of the FGM by evaluating time-ordered exponential 
functions. Equivalent but more complicated forms of this equation were first derived 



where i,j = 1, 2 corresponds to 



^In this work we will identify matrix elements as Uij with Ua 
a,a' = +,-, e.g. U12 = U+-. 

*Note, that this definition deviates from the definition used in Ref. |24| . Here, the vectors and 
Va are given by the first and second column of the matrix Ua multiplied from the left by ±0-3 . 
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by Abrikosov and Ryzhkin Of special interest will be the trace of the Green 

function at coinciding space points, 

^n^G {x,x;uj)]^ — , , , . ■ (85) 

U+22{x)Ul22{x) - U+i2{x)Uli2{x) 

It immediately follows from Eq. (^4|) that the local DOS is given by 

, X 1 f U+22{x)Ul22ix) + U+12ix)Ul,2{x) \ 

^ ^""^ W22ix)Ul,,i.)-U,MUl,2ix)) ■ ^''^ 

3.3.1 Riccati equation 

Since Eq. (^) only depends on the ratios^ 

<^a{x) = U^12{x)/Uo,22{x) , (87) 

we may also write the DOS as 

p{x, u;) = -Re , ^^)';)' • (88) 



TT \1 - $+(x)$l(x) 

More generally, it follows from Eq. (^) that the whole matrix Green function at 
coinciding space points can be written in terms of $+(x) and $l(x): 

iG^ix, x;uj) = ^ [g^{x + 0+, x; u) + g"{x, X + 0+; lu)] 



\{l + ^+{x)^*_{x)) -^+{x) 



1 -$+(a;)$l(x) V -$l(a;) \ {l + ^+{x)^*_{x)) 



(89) 



Using Eq. (77), we find that the $c<(x) are both solutions of the same Riccati equation, 

- id^^a{x) = 2Co{x)^a{x) + A(x) + lS.*{x)^l{x) . (90) 

Here, we have introduced ij(x) ~ to — ^ix). Similar Riccati equations have recently 
been obtained by Schopohl [ p8[ from the Eilenberger equations of superconductiv- 
ity. To specify the initial conditions, let us assume that outside the interval [0, L] 
the potentials V{x) and A(x) are real constants, Vbc and Abc > 0. This amounts 
to taking the limit A ^ oo, but keeping L constant. The initial values <&+(0) = 
limA^oo 5i2(0, -A)/522(0, -A) and $_(L) = limA_oo(5-^)i2(i + A, L)I{S~^)22{L + 
A, V) can be obtained by evaluating the S'-matrix for constant potentials. 

3.3.2 The S-matrix for constant potentials 

For constant potentials Vn and A„, the time-ordering operator T may be omitted in 
Eq. (|8l|), and the 5- matrix is given by 



S„{X - X') = COSh[V|A„|2_^2 (.^ _ 2.')] 



+ zsinh[V|A„P-^g(.-x')] ■ (91) 



:,2 



^Note also that the definition of ^a{x) does not involve the extra factor ±i used in Ref. |24 
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For |A„p < <D^, the argument of the square root is negative, so that in this case we 
write the S'-matrix as 



Sn(x - x') = cos[v/iI'^ - |A„|2 (a; _ x')]ao 

+ I sin[VcD2-|A„|2 (x - x')] i^"'"' 



A„a+- A:,ct_ 



^^l - I A. 



(92) 



For notational simphcity, let us introduce A^^'''* = ^|A„p — cj^. To calculate $+(0) 
and $_(L), we take the limit x — > ±cxi of the ratio Sni2{x) / Sn22{x) and obtain 



lim 



Snl2{x) 



x^±oo Sn22{x) 



(93) 



If we keep in mind that the frequency uj involves a small imaginary part, we see that 
this result is not restricted to the case |A„p > ci^: For |A„p < the square root 
has to be taken such that the right-hand side of Eq. ( p3| ) vanishes as A„ 0. This 
follows directly from the definition of the S'-matrix. 

3.3.3 Initial conditions 

It follows from Eq. (|9^) that the Riccati equation (|9^) should be integrated with the 
initial conditions 



$+(0) = $1(L) 



Abc 



sgn(^-yBc)V(^~VBc)^-A|c-(^-"t43c) 



A|c > - Vbc? , 
A|c < ('^ - ^BC)' , 



(94) 



where sgn(x) is equal to 1 for positive x and equal to —1 for negative x. While for 
Abc = the initial conditions are given by $+(0) — ^*_{L) = 0, for Abc ^ oo one 
gets $+(0) = $1(L) = i. Note that for arbitrary potentials Vbc and Abc: the initial 
values are simply given by the stable stationary solution of the Riccati equation ( |90| ) 
with V{x) = Vbc and A(a;) = Abc- 



3.3.4 The case of a discrete spectrum 
Defining the (complex) phase (pa{x) via 

$„(a;)=e*'^-(^) , 

and decomposing A{x) into its amplitude |A(a;)| and its phase {}{x) 

A{x) = |A(a;)|e*''(") , 
the Riccati equation (pG) turns into 



dxfaix) = 2u){x) + 2|A(a;)| cos [(pa{x) — ^{x)] 



Note that this equation of motion is of the Langevin type 59, 60[ 
dxv{x) = -a + biV{x) + b2{v)A{x) + 52(w)A* {x) . 



(95) 
(96) 

(97) 

(98) 
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Let us consider the case (w — Vbc)^ < ^bc- directly follows from Eq. (Q) that 
I $+(0)1 = |<i>_(L)| = 1, such that the initial values (p+{0) and f-{L) are real. Hence, 
the solutions of Eq. ( |97| ) remain real, which implies |$q,(x)| = 1 for all x. As can be 
seen from Eq. (^), and ^P-{L) can be chosen to fulfill '^+(0) = —(p-{L) € [0, tt], 

so that the initial values V'+(0) and ip-{L) are uniquely determined by 



cot ip^ (0) = — cot (L) = 



BC 



(w - VBcy 



(99) 



For the phases (paix) to be continuous, they have to be unreduced phases which are 
not limited to take values between and 2tt. In terms of the (pa{x), the local DOS 
can be written as 



1 

-— Imcot 

TT 



2 6 {(p+{x) — (p-{x) — 27rm) 



(100) 



m— — oo 



It is easy to show that for x > the phase (p+{x, lu) is a monotonic increasing function 
of uj. For the initial value at a; = 0, this follows from Eq. (99): 



1 



BC 



> 



(101) 



It can be seen directly from Eq. (|97|) that d^ip+(x) > is true for every x > 0. More 
formally, differentiating the equation of motion (|97|) with respect to uj gives 



d.j,d^ipa{x,uj) = 2 - 2|A(a;)| sin[(/7Q(a;) - ^{x)] d^ipa{x,uj) . 
Solving this first-order differential equation for dt^(pa{x,uj), we obtain 

d^iPa{x,Uj) = d^ipa{xo,Uj) 

\A{x")\sm[ipa,{x")-i9{x")] dx" 



(102) 







/ exp 




Ixo 





dx' 



(103) 



Since both terms on the right-hand side of Eq. (103) are positive for a = + and 
Xo = 0, we have d^'pj^{x,u}) > for a; > 0. Analogously we find dt^ip-{x,u}) < for 
X < L. For arbitrary uj, ip-^-{x,uJ) integrated with the initial condition <y9+(0,a;) given 
in Eq. ( p^ ) will usually not be equal to (p-{L,uj) up to a multiple of 27r at a; = L. 
For certain discrete frequencies ujm, however, this is the case. Since (/3+(a;,w) is a 
monotonic function of w, we can uniquely define ujm by the condition 



ip+{L,uj,-n) = ip-{L,ujjn) + 27rm . 



(104) 



Note that ujm is only well-defined if it turns out that {uJm — Vbc)^ < Since the 

right-hand side of Eq. ( p7[ ) is a 27r-periodic function of ipa{x), we see that ip+{x,ujm) 
and f-{x, UJm) are equal up to the constant 27Tm for every x € [0, L], 



ip+{x,ujm) - Lp-{x,ujm) = 27rm 



(105) 
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Of course, the ujrn are the discrete eigenvahies of the system. This follows immediately 
from the fact that the local DOS p{x, lu) is equal to zero if not w = ujm for one to. 
Therefore, the total DOS is given by 



(106) 



Using the well-known formula 6{f{uj)) = \f'{uj, )| '^(^ ^ ^m), where ujm are the 
zeros of f{oj), we can write Eq. (100) as 



2S{UJ - UJrn) 



(107) 



Expressing dujfa{x,uj) by the right-hand side of Eq. ( |103D with a = -I- and xo = or 
a = — and xq = L, respectively, we arrive for Abc = oo at 



dx' exp (-2 dx" I A(x")| sin [(^„(x", c^™) - i3{x")]) 



(108) 



Once the eigenvalues LUm have been determined, this equation in principle allows to 
calculate the local DOS for arbitrary potentials V{x) and A{x). 

3.3.5 Integrated averaged Green function T{lu) 

In the last subsection, we have seen that the integrated DOS can be obtained by 
solving a simple initial value problem for the phase 1^3(2;), which is a functional of the 
disorder. To implement the correct boundary conditions for a system of length L, we 
have first assumed that outside the interval (0, L) the potentials are constant over a 
range A, but then drop to zero. Finally, we have let A go to infinity. However, if we 
are only interested in the limit i — > 00, i.e. the bulk properties, we can set A = at 
the beginning of our calculations. The physical meaning of this is that bulk properties 
should be independent of the boundary conditions. 

By setting the potential equal to zero outside the interval (0, L), we will not only 
be able to recover the equation of motion (97) satisfied by the phase (p{x) which 
determines the integrated DOS, we will also be able to derive an additional equation 
which allows to calculate the inverse localization length. 

It follows from Eq. (|8^) that the trace of the space-averaged diagonal element of 
the retarded Green function is given by 

(Tr[g«(x,x;u;)])^ = j ["^ dx Tr[g«(x, x; a;)] 
1 



iS22(i,0) 




L 

i I dx {S22iL,x)S22{x,0) - S2l{L,x)Si2{x,0)) 




(109) 



The term in angular brackets can easily be identified to be equal to d^S22{L, 0). We 
can therefore rewrite Eq. (109) as 



(Tr[c;«(x,x;a;)])^ 



d^S22{L,0) 
LS22{L,0) 



^a^ln[522(i,0)] 



(110) 
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To describe the bulk properties, one should now take the limit L — > oo. If we introduce 
r(u;) as the trace of the energy-integrated space-averaged Green function at coinciding 
space points, 



T{uu) = lim \\n[S22{L,Q;uo)] 



(111) 



(112) 



we have 

(Tr[5«(x,x;c^)]) =9^r(w) . 

While the integrated DOS is given by J\f{Lu) = — 7r^^Inir(a;), we will see in the next 
subsection that Rer(cL') is equal to the inverse localization length ^^"^(1^). The decom- 
position of T{uj) into its real and imaginary part can therefore be written as 

T{lu) ^ e-^{Lu)~inAf{uj) . (113) 

Both A/'(iu;) and £~^{u>) can simultaneously be calculated by determining the logarithm 
of the S-matrix element S22 ■ It is convenient to express the S'-matrix elements in terms 
of their phases. Let us define ipaa' (x) via 

5„„.(x,0) = e-*'^"='(^) . (114) 

r(a;) is then given by 



iT(uj) ~ lim 

L- 



(115) 



Introducing a = —a, the properties of the S'-matrix Saa' — S*^^, render into ipaa' = 
— The S'-matrix can be expressed in terms of ipi2, ^22 and its complex conju- 
gates.H It follows from idxS{x,0) — A-'I{x)S{x,0) that the ipaa' satisfy 



dxfaa'ix) = Maaix) + M^aix) ejq) [i(<y5QQ'(x) - (faa'ix))] 



(116) 



Recalling that M(x) 
(f22 read 



-'Lu(x)aT, — A(a;)(T-|_ + A*(a;)(T_, the two equations for (pi2 and 



dx<P22ix) = Ujix) + A*{x) exp [i{ip22ix) - if\2[x))\ , 

dx^\2{x) = -u}{x) - A(a;) exp \-i{ip22{x) - ^\2{x))\ 
If we now introduce 

ip{x) = (p22ix) - ipi2ix) , 
C(x) EE ~i {ip22{x) + (Pl2{x)) , 

we arrive at the following system of equations of motion: 



dxVix) 


= 2u}{x) - 


f 2|A(a;)|cos[(^(a;) - i9(a;)] , 


dxC{x) 


= 2|A(x) 


sin [1^9(2;) — i9(x)] . 



(117) 
(118) 

(119) 
(120) 



(121) 
(122) 



"Recall that instead of a, a' = -|-, — we also use i,j = 1,2 (see footnote M on page |20[) 
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Note that Eq. (121), which determines ip{x), is exactly the same Langevin equation 
that we derived before from the Riccati equation [see Eq. (|97|)] and is independent 
from Eq. (122). After having found a solution to Eq. (121), ((x) can in principle be 
obtained by integrating Eq. (|122|). In terms of (p{x) and C{^)> r(aj) is now given by 



iT^Lu) ~ lim 

L — *oc 



ML) + tC{L)]/2L 



(123) 



The initial condition S{0, 0) = ctq could be mapped on the initial conditions for (p{0) 
and C(0) which are, strictly speaking, singular but integrable. In the limit L —^ oo, 
the initial conditions finally drop out, but to describe bulk properties for finite L, it 
is better to choose ip{0) = C(0) = 0, such that Lp{x) and ({x) are real for all x and 
there is no contribution to r(w) which as L becomes large only dies out as 1/L. The 
integrated DOS and the inverse localization length can now be expressed as 





Po lim Lp{L)/2L , 




= lim CiL)/2L . 







(124) 
(125) 



These two equations in combination with the equations of motion (121) and (122) allow 
for simultaneous exact numerical computations of the (integrated) DOS and the inverse 
localization length for arbitrary given disorder potentials. Since the (integrated) DOS 
and the inverse localization length are self-averaging quantities it is sufficient 
to consider just one typical realization of the disorder potential. We will do this for 
various interesting cases in Section [sj 

In analytical calculations one does not usually work with a certain realization of 
the disorder. Instead, one tries to calculate averaged quantities by using the given 
statistical properties of the disorder potentials. Taking the average of Eq. (123) with 
respect to the distribution of the random potentials V{x) and A(x), we obtain 



i{T{Lu)) — lim 

L — >cc 



mL))+t{C{L))]/2L 



(126) 



Integrating the equations of motion (121) and (122) with respect to x from to L, 
(r(cij)) can be rewritten as 



i{T(uj)) = lim 



1 

OO L 



dx {Cj{x) -f A* (a;) exp . 



(127) 



The average in Eq. (127) has to be taken with respect to the probability distribution 
involving the disorder at x and via (p{x) also at all space points between and x, 
which can also be considered as an average with respect to the joint probability dis- 
tribution of the random potentials at x and the unreduced ^hase ip{x). However, since 
exp [«v'(2:)] is a 27r-periodic function in ip{x), it is also sufficient to use the joint proba- 
bility distribution of the random potentials and the reduced phase ip{x) e [0, 27r). For 
large x, the joint probability distribution becomes stationary, and since {V{x)) — 0, 
Eq. (p7|) reduces to 



i (r(cj)) =uj+ {A*{x) exp [iip{x)]) 



(128) 
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We will use Eq. (|T2|) in Section | to calculate both M{uj) and ^^H^) for the FGM in 
the white noise limit. 

3.4 Gauge invariance 

It turns out that fluctuations of the forward scattering disorder have similar effects 
on the DOS and localization length as have phase fluctuations of the gap parameter. 
To illuminate the hidden symmetry, let us again consider the equation satisfied by the 
retarded Green function {x,x' ;u!), 

A crucial point of this equation is that its form is left invariant under the gauge 
transformation pl| 

x'; Lj) ^ exp [-{t/2)x{x)a3] g^{x, x';lo) exp [(z/2)x(a;')^3] , (130) 
Vix) ^ Vix) + '^d^xix) , (131) 
A(a;) ^ A(x)exp[-ix(a;)] . (132) 

Here, x(a;) is a local phase rotation which is allowed to vary arbitrarily from point to 
point. Since the trace of a product of matrices is invariant under cyclic permutations 
of the matrices, Tr g^{x, x; uj) and therefore the local DOS and the inverse localization 
length are also invariants under the above gauge transformation.]^ 

Instead of directly looking at the trace of the Green function at coinciding space 



points, we can also consider the equations of motion (121) and (122). Their form is 
gauge invariant under the combined transformation 

ip{x) Lp{x) ~ x{x) , (133) 

C(x) ^ Cix) , (134) 

Vix) ^ Vix) + ^d^x) , (135) 

i^ix) dix) - xix) . (136) 



It follows from Eqs. (124) and (125) that both the integrated DOS A/'(w) and the inverse 
localization length i^^iuj) are invariant under the considered gauge transformations. 
Only in the unusual case of a finite limit lim^^oo xi^)/ L we get an irrelevant shift in 
the additive constant of A/'(w). The choices for which either Vix) or dix) vanishes are 
especially convenient. 

1. Effectively vanishing phase fluctuations: If i?(x) is differentiable, we can define 
xix) = "dix), such that there are no phase fluctuations of A(a;) left and A(a;) 
can be taken to be real and positiv. We will use the resulting transformations 
Vix) Vix) + dx'dix)/2 and A(x) |A(x)| at the end of Section | to find an 
exact solution for the FGM involving only phase fluctuations. 

^Strictly speaking, it follows only that d^l~^{u}) is left invariant. The integration constants at 
a; = CO are, however, the same, so that l~^{u}) is invariant under the gauge transformation. 
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2. Effectively vanishing forward scattering potential: Choosing the phase xi^) such 
that dxX{x)/2 = —V{x), the forward scattering potential V{x) can be eliminated 
by renormalizing the phase fluctuations 'd{x) of the backscattering potential 
A(x). 



3.5 Lyapunov exponent and localization length 

In this subsection, we will explicitly show that the Thouless formula holds for the FGM 
and that Rer(cL;) can indeed be identified with the inverse localization length £^^{uj). 

One of the striking properties of disordered systems treated in the independent electron 
approximation is the fact that the disorder can lead to a macroscopic large number 
of localized states. These states are eigenfunction of the Schrodinger equation falling 
off exponentially with distance from one point in space which is characteristic for the 
particular solution. In one dimension, an arbitrary weak disorder suffices to localize 
all eigenstates (excluding perhaps states at isolated energy values) |6|]. As a conse- 
quence, the diffusion coefflcient and the dc conductivity vanish. Therefore, there is no 
metal-insulator transition in one dimension. 



3.5.1 Thouless formula 

Since the energy dispersion of the FGM is linear, the Schrodinger equation of the FGM, 
Hip(x) = ujip{x), is a linear first order differential equation. Fixing the two-component 
wave function i^i^) = {'4'i{x),ip2{x))'^ at one space point xq therefore constitutes the 
wave function at all space points x. As we will see explicitly below, for large distances 
— a;o| the envelope of the wave function will grow exponentially with probability one, 
i.e. ||V'(a;)|| ~ ||V'o|| exp (-H7|x - xqI), where ||V'(2^)IP = \^iix)\'^ + \^2{x)\'^ ■ Of course, 
ipix) cannot be an eigenfunction of the Hamiltonian H satisfying the right boundary 
conditions. The proportionality factor 7 in the exponential function is called the 
Lyapunov exponent. The (mean) localization length is usually defined as the inverse 



Lyapunov exponent |20|, 



rV)..H..»J..(M£f). (137, 

Note that since we take the limit L ^ 00, the initial value \\tJJo\ \ drops out. 

Let us now rewrite the Schrodinger equation such that the differential operator dx 
is proportional to the unit matrix. As in the non-Abelian Schwinger-ansatz we factor 
out a (73-matrix so that the wave function ipi^) = o's'^ix) satisfies 

[idx - M{x)] ip{x) = . (138) 

Here, M{x) is the matrix defined in Eq. (^3). The solution to the Schrodinger equation 



(138) is therefore given by 

?/i(x) = S{x,xo)ipo ■ (139) 

It follows with ||'0(a;)|| = [[-(/"(a;)!! that the Lyapunov exponent can be expressed in 
terms of the S'-matrix, 

7(cj) = lim i In S{L, 0) i^o ■ (140) 

L — ^OQ _L 
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Due to l^iip = |5i2|2 - 
equally fast and we have 



1, S22 = 'S'li and 5*21 = S12^ all 5- matrix elements grow 



1 



7H= lim -ln|522(i,0)| 

L— »oo ij 



(141) 



Comparing this equation with Eq. (Ill), we see that the inverse localization length 



£ ^(w) = 7(ti;) is in fact equal to Rer((jj). The equation 



r^{Lo) = Rer(w) 



(142) 



is known as the Thouless formula pd^ and was first shown to be valid for the FGM 
by Hayn and John |]67| in a different and more complicated way. But note that these 
authors have only derived Eq. (142) up to an integration constant which they had to 
determine by different means. 



3.5.2 Localization length at uj — for real A(x) 

Although we are going to postpone detailed calculations of the DOS and the localiza- 
tion length to the next sections, let us now consider the localization length at frequency 
w = for a real disorder potential A(x) and V{x) = 0. In this case, M{x) — —ia2A{x), 
so that the S'-matrix may be expressed without the path-ordering operator as 



S{x, x') = cosh 



A(y) dy 



(To ^ sinh 



A(y) dy 



<J2 



(143) 



It follows from Eq. (141) that the inverse localization length at cj = is given by 



"1(0) = lim -In cosh 



A(x) dx 



= |(A(x)),| = |A, 



(144) 



where we have used the fact that the average (A (a;)) a: = lim/,_>oo 7; /q^ ^(2:) is equal 
to the expectation value Aav The inverse localization length at frequency w = is 
equal to the absolute value of the expectation value of the backscattering potential 
(A(a:)) and does not depend on the random fluctuations around this average value. 
Note that this result is valid for arbitrary higher correlation functions of A(a;). In the 
case (A(x)) =0, the localization length diverges for w = which clearly distinguishes 
the point w = 0. The uj — Q eigenstate is delocalized! 



4 Exact results 

In this section, we review exact results for the density of states (DOS) and the inverse 
localization length of the fluctuating gap model in the white noise limit and the limit 
of infinite correlation lengths ^. 

*It should be noted that the above reasoning is not true in the unprobable case where i/>o is very 
close to the eigenvector to the smallest eigenvalue of S. However, since the definition of the Lyapunov 
exponent involves the limit L — > 00, this only happens with zero probability. 
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4. 1 The white noise hmit 

For small correlation lengths ^, the disorder of the fluctuating gap model (FGM) may 
be approximated by Gaussian white noise. This ^ — > limit is of special interest 
because in this case the disorder at different space points is uncorrelated which basi- 
cally admits for an exact analytic solution of the model. Various methods may now 
be applied to find analytic results for the density of states (DOS). Ovchinnikov and 
Erikhman [p^ were the first to solve the commensurate case for which the random 
backscattering potential A(x) is real. They showed that in the symmetric phase for 
which (A(a;)) = 0, the DOS has a Dyson singularity previously found by Dyson 
In the case of (A(x)) ^ which models a phase below a phase transition, the DOS 
either exhibits a singularity or a pseudogap near the Fermi energy depending on the 
ratio of the disorder and the static gap. Using the technique of S'-matrix summation, 
Golub and Chumakov |Q| confirmed the results by Ovchinnikov and Erikhman and 
were also able to solve the incommensurate case. In the incommensurate case which 
is described by a complex backscattering potential, there is no singularity and the 
disorder can only lead to a filling up of the pseudogap. The incommensurate case was 
also considered by Abrikosov and Dorotheyev . 

In recent years, the method of supersymmetry developed by Efetov [ |66[ has been 
established as a powerful tool to describe disordered systems in the white noise limit. 
First, Hayn and John JSTf re-derived the Ovchinnikov and Erikhman result for the 
DOS and were also able to give an analytic expression for the localization length. 



Later, Hayn and Fischbeck [^,69| used the method of supersymmetry to generalize 
the known results for the integrated DOS to the case of three independent disorder 
parameters which describe forward, backward, and umklapp scattering. These solu- 
tions include both the commensurate and the incommensurate case as special cases. 
Finally, both the integrated DOS and the localization length were calculated by Hayn 
and Mertsching in the most general case with a complex static gap parameter and 
three disorder parameters. 

In this section, we use the formalism developed in the previous section and follow 



the ideas of Lifshits, Gredeskul and Pastur |20| to show that the probability density for 



the distribution of the reduced phase ip satisfies a continuity equation. The stationary 
probability flux of the continuity equation turns out to be equal to the integrated DOS. 
We then derive an equation closely related to a Fokker-Planck equation which allows 
to rederive the integrated DOS for the most general case exactly. Before considering 
the most general case we discuss the commensurate case, i.e. the Ovchinnikov and 
Erikhman limit which we can solve in analogy to Halperin's calculation of the inte- 
grated DOS of a particle with an effective mass in a white noise disorder potential [ [70| . 
The treatment of the general case is similar but more awkward than the Ovchinnikov 
and Erikhman limit because instead of a linear differential equation of second order 
one has to face a linear differential equation of fourth order. The equations to deter- 
mine the integrated DOS and the localization length are, however, the same as those 
derived by Hayn and Mertsching using the method of supersymmetry |^ so that we 
recover their general results which also include the incommensurate case which we will 
discuss afterwards. 
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4.1.1 Equality of the integrated density of states and the stationary probability Rux 
As shown in the previous section, the averaged integrated DOS can be written in the 
thermodynamic hmit as [see Eq. (|128|)] 



U{u;)^{FUV,A,^))/2tt . (145) 
F^{V, A, (f) is hnear in the disorder V, Re A, Im A and given by 

F^(y,A,ip) = 2(cj- V) + 2ReAcos^ + 2ImAsin(^ . (146) 



The unreduced phase ip{x,uj) satisfies the equation of motion (121) which can also be 
written as 



dx(p{x,uj) ^ F^{V{x),A{x),ip{x,uj)) 



(147) 



Continuity equation 

The space-dependent probability distribution can be defined as 

Puj{x, ip) = {52T,{(p - (p{x, LU))) , 



(148) 



where S2ti{x) = X]m=-oo'^('^ ~ 27rm) is the 27r-periodic delta function. A continuity 
equation may be derived by partially differentiating P^^ix, tp) with respect to x which 
results in dxPuj{x, p) — -~d ^{52Tr{p — ip{x,ijj))dx'p{x,uj)) . Replacing dxP^{x,Ljj) by the 
right-hand side of Eq. (p_47[) , the continuity equation reads 



(149) 



dxPwix, ip) + d^J^{x, if) ^ 



where the probability flux Juj{x, ip) is given by 

Juj{x,(p) = {S2Trip - (pix,uj))F^{V{x),A{x),ip{x,uj))) 



(150) 



Letting x go to infinity, the probability distribution and the probability flux become 
stationary, i.e. independent of a;. Due to the continuity equation (149), J^^ becomes 
also independent of p. Integrating the stationary form of Eq. (15C) with respect to (p 
from to 2tt therefore leads to 2ttJ^ — {Ft^{V, A, p)), so that together with Eq. (145) 
we find the remarkable relationship EO] 



(151) 



The integrated DOS, i.e. the number of states in the energy interval (0,Ci;] per unit 
length, is equal to the stationary probability flux. Note that this result is valid for any 
disorder potential and is not restricted to the white noise limit which we will consider 
in the following. 
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4.1.2 White noise and the Fokker-Planck equation 

While for arbitrary finite correlation lengths ^ and Gaussian statistics, it does not 
seem to be possible to find an exact analytic expression for the (integrated) DOS, the 
white noise limit ^ ^ 0, V^^ Dy, (ReA^)^^ Dr, and (ImA,^)^^ Dj admits 
for an exact solution. This is due to the fact that in this case the disorder at different 
space points is uncorrelated. In the white noise limit, the disorder is characterized by 
the correlation functions 

{V{x))=0, {Vix)V{x')) =2DvS{x-x') , (152) 
(ReA(a;)) = ReAo , {Re Aix)Re A{x')) = 2Dii6ix - x') , (153) 
(ImA(a;)) = ImAo , {Im A{x)lm A{x')) 2Di S{x - x') , (154) 

where A(a::) = A{x) — Aq. It is no loss of generality to assume that the first moment of 
the forward scattering potential V{x) vanishes because a finite value would only lead 
to a renormalization of lu. Since the probability distribution of the disorder is assumed 
to be Gaussian, higher correlation functions arc simply given by Wick's theorem. 

To cast the continuity equation into a Fokkcr-Planck equation, we make use of the 
Gaussian nature of the disorder, so that for a functional f{V(y)} of the disorder V{y) 

we have {V (x) f {V (y)}) = J dx' {V{x)V{x')) {^^4F(iPy")' "^ti^rc in the last term we 



take the functional derivative of f{V{y)} with respect to V{x'). Using Eq. (152), this 
simplifies to 

{V{x)f{V{y)}) = 2Dy 1^-^^^) ■ (155) 



If we want to apply this relation to Eq. ( |150| ), f{V{y)} has to be replaced by the 
27r-periodic delta- function 52'r^{f — (p{x,uj)) whose phase ip{x,uj) is a functional of the 
disorder involving the disorder at all space points y < x. Using the chain rule for the 
functional derivative, we get 



{S^Avix, OJ) - V{x)) = -2Dvd^ ( S^A^ix, oj) ~ ^) ) . (156) 



Writing '.p{x,Lu) as (/?(a;,u;) — (/3(0,w) = Jq dx' F^{V{x'), A{x'),ip{x' ,uj)), we find 
^svlx^) ~ ~26'(a; — x'). Since in our case x — x' , 9{0) needs to be defined care- 
fully (see also Itzykson and Drouffe |7l[]). Recalling that we have introduced the 
delta function S{x) as the limit ^ — > of a symmetric function of x, we see that 
to maintain this symmetry, we have to define 9(0) = 1/2. It therefore follows with 
Puj{x,ip) = {S2-^i^p - (p{x,LLi))) that 

(<52x {^~^{x,uj)) V{x)) - 2Dvd^P^ (x, if) . (157) 

Similarly, we can show that 

{d2Tr{^ ~ ^ix,u!))Re A{x)) = -21)^9^ (cos </? P^, (x, (^s)) , (158) 
(S2Av-vix,uj))lTnA{x)) ^ -2DidAsmipP^{x,ip)) . (159) 
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Making use of these relations, the probability flux Ju>{x, ip) is given by 

Juj{x, ip) = 2[uj + Re Aq cosip + Im Aq si-aLp\P^{x, Lp) 

- 4 \Dvd^P^{x,ip) + Djicosifid^ {cos (pP^{x, Lp)) 

+ Di sin ipd^ {sin tpPuiiXjip))] . (160) 

The probability flux can also be written as 

JUx,ip) = A^{ip)PUx,ip)-^d^[B{ip)PUx,ip)] , (161) 

where 

Auj{'p) = 2[uj + Ke Aq cos ip + Im Aq sin ip] — 4(1?^ — Dj) cos (p sin ip , (162) 
B{(p) ^8[Dv + DRCOs^ip + Disin^ip] . (163) 

Note that the forward scattering disorder Dy only leads to a renormalization of Dn 
and Dj. We therefore define Da = + Dy and Dj = Dj + Dy- 



Eq. (161) together with the continuity equation explicitly shows that the probabil- 



ity distribution satisfies the following one-dimensional Fokker-Planck equation |59|, |60[ 
dM^,^) = [AUv)PU^,v)] + [i?(^)P^(x,^)] . (164) 



To find an analytic expression for the integrated DOS, we use the fact that, as shown 
above, M{uj) is equal to the stationary probability flux. A good starting point to 



calculate the integrated DOS is therefore the stationary form of Eq. (161), which is 
nothing but the integrated stationary Fokker-Planck equation with the constant of 
integration being the stationary probability flux that is equal to the integrated DOS 

M{io) = A^{^)P^{^) - [B{^)P^{^)] . (165) 

In principle, one could flrst flnd a solution for [5(1^9)^1^(1^9)] to this first-order differ- 
ential equation subject to the boundary condition [_B(27r) P^{2tt)] — [B{0) Puj{0)] with 
/^{uj) as a parameter and could then use the normalization condition ^ dip Pt^{ip) — 1 
to determine Af{uj). While this procedure works quite well in the incommensurate 
case with Dji — Dj and after a variable transformation also in the commensurate case 
(without forward scattering) for which Dj = Dy = 0, the most general case consid- 
ered here seems to defy such a treatment. Here, we therefore present an alternative 
method which allows to recover all known results in the white noise limit. 

4.1.3 The density of states in the white noise Umit 

To find the (integrated) DOS for arbitrary parameters Djf ,Dj,Dy and complex Aq, 
we first use the variable transformation 



z — tan 



(l-f). («e) 



The trigonometric functions sin tp and cos tp can now be expressed in terms of z as 
i^i^ , and cos (p — ~ ■ Since 4^ 

1+2^ ' r l-\-z^ dtp 



simp = , and cosip = ~t^- Since = 5(1 + z^), we find = i(l + z^)dz 
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and the probabihty distribution P{(p) has to be replaced by i(l + z^)P{z). Using 
these relations, Eq. (165) turns into 



Miuj) = (cj + ImAo) - 2(ReAo - {^Dr - Dj))z + {uj - IniAo)^^ 



P{z) - bid:, ( 1 + 2DY\2Db, - bi)z^ + P{z) 

Taking the Fourier transform of this equation leads to 

27r7V(w)(5(fc) = (w + Im Ao)P(fc) - 2i(Re Ao - (2I)h 
- (w - Im Ao)P"(fc) - 2ibiP"'{k) 



(167) 



iDik 



P{k) - 2DY'^{2Dr - Di)P"{k) 



Dj))P'{k) 
P""{k) , 



where 



P{k) = / e-''^^P(z)dz 



(168) 



(169) 



is the Fourier transform of P{z) which is also known as the characteristic function [ [72[ . 
Normalization of the probability distribution P{z) implies P{k = 0) = 1. Assuming 
higher derivatives of P{z) to be integrable, the other boundary conditions are given 
by demanding that P{k) — > sufficiently rapidly as |fc| — > oo. 

4.1.4 The commensurate case without forward scattering 

Before we proceed with the most general case, let us first consider the Ovchinnikov 
and Erikhman limit, i.e. the commensurate case without forward scattering. In this 
case, Aq is real and may be assumed to be positive, D/ = Dy = and we set D = Dr. 
Eq. (168) reduces to a second order differential equation with the boundary conditions 
P(0) = 1 and P{k) ^ as |fc| ~* oo: 



2TTAf{Lu)6{k) = LjPik) + AiD 



2D 



Integrating this equation from — e to +e with e 



27r7V(w) -tj P'(0+) - P'(0") 



AiD 
0+ leads to 



Since P{z) is real, P{—k) — P*{k), and P'{—k) = -P'*{k), which implies 



ReP'(0^ 



(170) 



(171) 



(172) 



P'(0+) can be determined by first finding a solution to the differential equation ( |170| ) 
for fc > which vanishes as k approaches infinity and then normalizing this solution 
such that P(0) = 1. Because Eq. (|17C| ) is homogeneous for fc > 0, it therefore follows 
that if g{k) is any solution to the homogeneous differential equation which obeys 
g{k) ^ as fc ^ oo, then 7V(tj) = Re [g' {0)/ g{0)]. We introduce y{t) = g{k) with 
t — —i{k + iut/AD). The integrated DOS is now given by 
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where y{t) has to satisfy the differential equation 



AD 



y{t) = 



(174) 



with the only restriction that y{t) should approach zero as t goes to ioo. The general 
solution of this differential equation can be expressed in terms of a linear combination 
of Bessel functions of the first and second kind, J,j{x) and N^{x) (Abramowitz and 
Segun (A&S) @). The only solution which satisfies the boundary condition involves 
the Hankel function of the first kind, H^\x) = J^{x) + iN^{x). Introducing 

Ao 



(175) 



we find y{t) = {vjlfu^j-^ [{oJt/D)^ 



/2 



Because the prefactor is real and we only 



need the imaginary part of the quotient y'[ijj/AD)/y[oj/AD), Eq. (173) turns into 



N{lo) = - Im 

TT 



(176) 



Note that this equation is valid for arbitrary lu and is even analytic in the upper half 
plane. In the following discussion of the Ovchinnikov and Erikhman limit, we will 
again restrict ourselves to u; > 0. It follows 



AA(w) 



•^^ {20) ^'v iw) -^i^ (20) (20) 



(177) 



The numerator can be simplified by using the Wronski relation [A&S, Eq. (9.1.16)] 
Ji, (x) Nl {x) — Jl {x) {x) ~ so that as our final expression for the integrated 
DOS we are left with 




(178) 



This result was first obtained by Ovchinnikov and Erikhman in a more complicated 
manner. Differentiating Eq. ( |178| ) with respect to uj, we find 



(179) 



where the derivatives of Bessel functions could in principle also be written in terms of 
Bessel functions. For Ag = which implies 1/ = 0, D is the only characteristic energy 
scale, and it is useful to measure all energies in terms of D. The insert at the top of 
Fig. shows the DOS p(lu) plotted versus uj/D. The DOS clearly exhibits a singularity 
near uj — whose asymptotic behavior can be found by noting that Jo(0) is finite while 
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Nq{x) diverges logarithmically for small x. It follow^ with Nq{x) ^ (2/7r) In a; [A&S, 
Eq. 9.1.8] 

such that the asymptotics of the DOS for Ao = is given by 

p(uj) ^ . (181) 

Recall that to generalize this result towards arbitrary frequencies w, we have to replace 
uj by The singularity described by Eq. ( |l8l|) is called a Dyson singularity and was 
first found by Dyson in a different model |0| involving also off-diagonal disorder. 
Outside this singularity, the DOS is almost equal to the DOS of the disorder-free 
model, taking its minima p{lo*) — 0.9636 po at lo* — ±1.25141?. 

If Aq 7^ 0, Aq is another characteristic energy scale, v = Ao/2D then basically 
gives the ratio of the two relevant energy scales Aq and D. The DOS plotted against 
w/Aq for different values of the parameter z/ is shown at the top of Fig. 0. While the 
singularity only survives for ly < 1/2, the DOS is constant for i/ = 1/2 and for v > 1/2 
the effects of the constant gap Aq dominate those due to the disorder and a pseudogap 
emerges. 

The algebraic dependence of the DOS at small uj on the parameter v can be found 
by using again an asymptotic expansion of the Bessel functions. If > is fixed and 

X ^ 0, the Bessel function J^{x) is finite and N^{x) {l/TT)r{i'){x/2)-'' [A&S, Eq. 

(9.1.9)]. It follows 

AD /■ UJ \'^^ 



SO that 



2iy / UJ \2i^- 



This implies that for i' < 1/2 the DOS in fact diverges algebraically as uj approaches 
zero because in this case the exponent is negative. Although the algebraic divergence 
differs from the divergence found in Dyson's model, we will nevertheless refer to this 
singularity as a Dyson singularity. For v > 1/2, however, the exponent is positive and 
the DOS vanishes algebraically. 

For large u the disorder becomes irrelevant and the DOS reduces to the mean-field 
result p{uj) = po e{uj^ - AI)uj/{uj^ ~ A2)i/2. 

^Note that the argument of the logarithm in the asymptotic form given by Ovchinnikov and 
Erikhman deviates by a factor 1/2 from our result. Nevertheless both expressions lead to the same 
asymptotic behavior. To take into account next to leading terms one has to use No{x) = 2/7r ln(aa;) + 
olx), where a = eT/2 which follows from A&S, Eq. (9.1.89). Here, 7 r; 0.5772 is Euler's constant 
which leads to a !^ 0.8905. This value is closer to our choice a = 1 than to Ovchinnikov's and 
Erikhman's choice a = 1/2. 
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Fig. 7 The DOS p(i^) and the inverse locahzation length I ^(oj) in the white 
noise limit for the commensurate case plotted versus oj/Aq for v — Ao/2D = 
0.1,0.3,0.5,1.0,3.0,10.0,100.0,00 or u = Ao/2D = 0.1,0.3,0.5,1.0,1.5,2.0,2.5,3.0, respec- 
tively. The inserts show the respecting graphs for the case Ao = 0, i.e. = 0. 
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Localization length 

Since r(w) = £~^(lu) — z7rA/'(a;) is an analytic function in the upper half plane, we can 
also easily find an analytic expression for the localization length £^^{u!). It follows 
from Eq. (176) that up to a constant 




The inverse localization length is now given by 



i~^iuj) = Rer(w) = -Lu- 



J2 J^] I 

^i' \2D) ^ V \2D / 



(184) 



(185) 



Comparing the right-hand side of this equation with Eqs. ( |178[ ) and (179), we find 

(186) 



LO p{lu) 



This equation is exact and can already be found Ref. [E0| 



If Aq = 0, it follows from Eqs. (180) and (181) that i ^{u) vanishes logarithmically 
as u approaches zero. 



2D 



In {uj/2D) 



Using Eqs. (182) and ( |183| ), we get for arbitrary Aq 
^-'(0) = Ao, 



(187) 



(188) 



which, as can be seen from Eq. (187) is also true for Aq = 0. Eq. ( |188| ) agrees with 
Eq. ( |144| ), so that T{ll!) involves no extra constant. For large frequencies, plco) pa 
and Af{uj) — > pa lu, such that 



D 



(189) 



A plot of the inverse localization length £ ^{uj) for various values oi v — Ao/2_D is 
given at the bottom of Fig. ^. 

4.1.5 Solving the general case with arbitrary parameters Dj^, Dj, Dy and Aq 
While in the commens urate case we only had to solve a differential equation of second 
order, for Dj ^ 0, Eq. (168) is a differential equation of fourth order and more difhcult 
to solve. Without loss of generality we may assume that Db. > Dj (later we can also 



take the limit Dji Dj). Integrating Eq. (168) from — e to +e with e O"*", we can 



proceed as before and express the integrated DOS in terms of P{k) and derivatives 
thereof evaluated at ±0^. Integrating the terms linear in k by parts and using again 
the fact that p(")(-fc) = (-1)"P(")* (fc), we find 



-Im 



-i{uj - Im Ao)P'(0+) + DjP"{0+) 



(190) 
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So, if y{k) is any solution to the homogeneous differential equation 

(to + ImAo)2/(fc) - 2i(Re Ao - {2Dr - Di))y'{k) - (w - Im Ao)2/"(fc) 

- 2ibiy"\k) - Wjk [y(k) - 2DJ\2Dr - Di)y"{k) + = , (191) 

which vanishes for fc — > oo sufficiently rapidly, then the integrated DOS is given by 



J\f{Lj) = -Im 



— — Im Ao) 



yio) ' y(o) 



(192) 



Again, since r(w) = £ ^{uj) 
up to a constant we have 



iTTj\f{uj) is an analytic function in the upper half plane. 



r(cj) = i{uj — Im Aq) 



y'io) 



(193) 



yio) y(o) 

Using the method of supersymmetry invented by Efetov ]6^ , Hayn and Mertsching 
p6| derived the set of equations (191) and ( |193D by different mean^^ Applying the 
method of Laplace transforms they found an exact expression with the constant of 
integration chosen such that one obtains the correct asymptotic behavior for large 
frequencies determined in the Born approximation, r(a;) = D — iuj = Dn + Di — icu. 
Here, instead of presenting a lengthy derivation, we will only cite the exact result 
found in I2II: 



T{uj) = 2Di + ADi 



F'(z) 

z{l - z) ^^^^ + zSr - i{l - z)e 



In this equation, F{z) is the hypergeometric hmction 

F{z) = F(i — ie + iSj — 6j^, ^ — ie — iSj — Sj^, 1 — 2ie; z) , 
with the parameters Sn, Si, e, and z (in our notation) given by 



ReAo 



AlDniDR-Di) 



1/2 



Si = 



ImAo 



1/2 



4 i^Di{DB, - Di) 

Dr-Di 
Dr 



1/2 



(194) 

(195) 
(196) 
(197) 



Recall that we incorporated the parameter of the forward scattering disorder, Dy, 
into Dr and Di by d efini ng Dr = Dr + Dy and Di = Dj + Dy . Only the additive 
constant 2Di in Eq. (194) does not get renormalized by Dy. This is due to the fact 
that for large frequencies T{lj) ^ Dr + Dj — ilo is independent of Dy . The imaginary 
part of Eq. ( |194| ) determines the integrated DOS, 



(198) 



^''Hayn and Mertsching use a slightly different notation but apart from this and some irrelevant 
different signs their expressions are equal to ours. 
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Taking the confluent limit Dj ^ of Eq. (194), one can recover Eq. (184) which 
describes the commensurate case without forward scattering. Turning on the forward 
scattering disorder Dy in the general expression gradually removes the singularity in 
the DOS for i/ < 1/2. Below, we are only going to discuss the incommensurate case. 
Other special cases can be found in 



4.1.6 The incommensurate case 

In the incommensurate case, Dr = Dj = D/2, and the forward scattering potential 
Dy only leads to a renormalization of D. If we intro duce D = D + Dy = Dr/2 + Dj /2 
and take the confluent limit Dj — > Dr in Eq. (194), we obtain 



r(cj) = D -iuj + Aa 



which is equivalent to 



-iuj/D 



(t) 



(199) 



T{lo)=D + Aq 



r 



iuj/D \ D 



-iu/D \ D 



(200) 



The imaginary part of Eq. (199) determines the integrated DOS, 




(201) 



This expression agrees with Ref. |Q. A plot of the DOS for different values of the 
parameter v = Ao/2D is shown at the left-hand side of Fig. |8[ There is no Dyson 
singularity, and in the absence of a static gap Aq, the disorder has no effect on the 
DOS so that p{u!) = po for — 0. At zero frequency, p{0) is always finite and the DOS 
vanishes with increasing as p(0) — po/ [/o(2i^)] . For a given Aq, the disorder leads 
to a filling of the gap. As in the commensurate case, in the limit Z? — s- 0, i.e. v oo, 
the DOS reduces to the mean-field result p(tj) ~ pa 9{uj^ — Aq) |cj|/(cj^ — A§)^/^. 
Taking the real part of Eq. (|19S|), we get for the inverse localization length 



r\uj) = D+Ao 



iui/D 



Ao\ T ~ (As>\ 
D J ^l-iuj/D \ D J 



T ~ Ao\ T - ( Ao 



(^) 



.(202) 



In contrast to the commensurate case, the localization length at w = is finite for any 
Aq and given by 



r\0) = D + Ai 



(203) 
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Fig. 8 The DOS p{io) and the inverse localization length £ ^{lo) in the white 
noise limit for the incommensurate case plotted versus u/Ao for u — Ao/2Z) — 
0.0, 0.3, 0.5, 1.0, 3.0, 10.0, 30.0, oo or i/ = Ao/2D = 0.0, 0.5, 1.0, 1.5, 2.0, 2.5, respectively. 



While for Aq/D — > one has ^^^(0) D (and also £^^{u!) D for every w), for 
Ao/i) > 1 one finds ^^^(0) ^ Aq + D/2 - Dv/2. For weak disorder, the disorder 
creates a few localized states with energies \uj\ < Aq whose inverse localization length 
is given by i~^{uj) ~ \/ Aq — u;^. A plot of the inverse localization length for different 
values of the parameter i/ and Dy = is given at the right-hand side of Fig. ||. 



4.2 InGnite correlation lengths 

In the limit of infinite correlation lengths, A(x) becomes independent of x and exact 
results for the averaged DOS of the FGM may be obtained. The limit of large cor- 
relation lengths ^ is of special importance for Peierls systems because the correlation 
length of the order parameter diverges at the Peierls transition. Sadovskii was the first 
to consider the fluctuating gap model (FGM) with infinite correlation lengths and 
calculated the one-electron Green function for the incommensurate case by summing 
up all diagrams in the perturbation expansion. This Green function leads to a DOS 
which exhibits a pseudogap at the Fermi energy. The commensurate case was later 
solved by Wonneberger and Lautenschlager ||7^ . 

Taking the ensemble average 

The limit of infinite correlation lengths can be solved by averaging the desired quantity 
calculated with a static gap A over an appropriate probability distribution of A. This 
amounts to taking an ensemble average. 



4.2.1 The commensurate case 

For real A and Gaussian statistics we have 
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Calculating the DOS by averaging 6{uj^ — A^) cj/V^^ — with respect to the above 
probability distribution, one obtains 



...M .^1^. /.(ijj. (205, 

Here, Iq{u) is the modified Bessel function with index 0. 

If we define the inverse localization length £^{oj) for ^ = cx3 by the Thouless 
formula, we have (^(uj) = (VA^ — uj'^ 6{A'^ — w^)), such that 

,2 poo 



^^'("V^^/ d-e-(-V2A^)"^y^:^. (206) 
4.2.2 The incommensurate case 

For complex A(x) and Gaussian statistics, the process of averaging can be written as 



dReA dImA 



(■■■)- J — — " e-|A|VA;._ (207) 
A similar calculation as above leads to 

(208) 



C'H=A.^e-""/i;. (209) 

Here, Erfi (u) = dx is the error function with an imaginary argument. 

Plots of Poc{^) and £l^{w) can be found in the next section. While in the com- 
mensurate case the DOS vanishes linearly in a;, it only vanishes quadratically in the 
incommensurate case. This is due to the fact that the probability distribution for com- 
plex A has less weight for small |A| than the one for real A. The inverse localization 
length assumes for both the commensurate and the incommensurate case a finite value 
at w = and drops to zero as uj increases. That ^^^(0) is finite in the commensurate 
case seems to contradict the general result ^^^(0) = Aav = derived in Section ||. 
One should keep in mind, however, that for ^ = oo we have only defined duji^{uj) 
by Ile{G{x,x;uj)). While for finite ^ a single chain is representative for an ensemble 
of chains, for ^ = oo, there is no self-averaging effect. On the other hand, it seems 
plausible to assume that the above results for ^ = oo give a good approximation to 
the case of finite ^ if f is much larger that any microscopic length scale involved. In 
particular, we have to demand A^^ ^ 1 and 3> 1. The above results for the DOS 
and the inverse localization length at w = can therefore not be expected to hold for 
finite correlation lengths. In fact, we will see in the next section that for any finite 
^ we find in the commensurate case p(0) = oo and £~^{0) = 0. For A^^ 3> 1 and 
uj^ 3> 1, however, we will find a remarkable agreement between the two solutions as 
predicted above. 
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5 Finite correlation lengths 



While in the limit of very small and infinite correlation lengths ^ of the random dis- 
order, the fluctuating gap model (FGM) admits for an exact analytic calculation of 
the density of states (DOS) and the inverse localization length, in the intermediate 
regime of finite ^ there are only approximate solutions available. It especially turns 
up the question: "How accurate are Sadovskii's solutions which for a long time 
were thought to be exact?" An answer to this question is of particular interest because 
Sadovskii 's solutions have become quite popular since the experimental discovery of a 
pseudogap in the underdoped cuprates above the critical temperature JI1|-[T|/. In 
this section, we will calculate the DOS and the inverse localization length for Gaussian 
statistics, as approximately done by Sadovskii with very high accuracy numerically. 
Finally, we will consider the case of only phase fluctuations for which we recently 
found an exact solution by applying a gauge transformation to the Green function and 
mapping the original problem onto a problem involving only white noise J2^ . 



5.1 Singularities in the density of states 

The exact results of the FGM derived in the white noise limit in the previous section 
imply under certain circumstances a Dyson singularity in the DOS. This singularity 
arises only in the commensurate case [i.e. for real A(x)] and only if the forward scat- 



tering potential and Aav = (A(a;)) are sufficiently small [see Eqs. ( |181| ) and ( |f83| )]. 
Since the white noise limit describes the low-energy physics of physical systems char- 
acterized by small correlation lengths this statement should also be true for small 
but finite f . As far as I know, it was first shown by myself in collaboration with Peter 
Kopietz that the DOS p{uj) of the FGM exhibits a singularity at the Fermi energy 
for any finite value of the correlation length ^ if the fluctuating order parameter field 



A(a::) is real and its average (A(a;)) is sufficiently small |22 . To detect the singularity, 
we applied the boundary condition Abc = Vbc = 0, such that the complete spectrum 
turned out to be continuous [see Eq. ( |9^ ) and its following remark]. In the Comment 
[77|, Millis and Monien showed that the local DOS p{llj = 0, x) calculated in our Letter 
[22| is equal to the absolute square of the wave function ip{x), which led them to claim 
that we have not calculated the DOS at all. However, as pointed out in the Reply 
[ [78| , for the boundary conditions used in Ref. [Q, one finds p{uj,x) = |V'w|^, i.e. the 
local DOS is equal to the absolute square of the wave function. This clearly invali- 
dates the argument given by Millis and Monien. Today, 1 would not use the artificial 
boundary conditions which lead to a continuous spectrum any more: The existence of 
the Dyson singularity put forward in the Letter [e3| can also be seen in the discrete 



case by considering the equation of motion (f21) for V{x) — and real A(a;) which 
after the shift Lp ^ f — 71/2 reads 

dxip{x) = 2uj + 2A(.t) sinLp{x) . (2f0) 

The Dyson singularity in the DOS is due to phase resonance: If lj is small (compared to 
As, A^^ and ^^^) but positive, the change of <p(x) is dominated by the fluctuating term 
2A(a;) sin 1^9(2;). Only near 1^9(2;) = wk (with n an integer) we have dx^{x) = 2lu > 0, 
such that (p{x) can only grow on average. As (<p(x) — tt-tt) w uj/Ag, fluctuation effects 
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of A(x) become important, driving 1^9(2;) from nvr + w/A^ to (n + l)7r — w/A^. Near 
(^(a;) = (n + l)7r, the constant force 2lj dominates again and the above picture repeats 
itself. 

As we decrease w, the "time" (which corresponds to the space coordinate x) to 
move 1^9(2:) from wk — oj/A., to titt + w/A^ will not change, but fluctuations of A(a;) 
will need slightly longer to drive ^p{x) from nir + w/Ag to (n + l)7r — to / /S.s, implying 
that 1^9(2:) decreases more slowly than w as w decreases. Now, the average DOS for 
frequencies between and lj is given by 

piCuj) = = hm , (211) 

where ^ is a number between and 1. Letting uj approach zero, it follows p(0) = 00. 
This divergence describes the Dyson singularity in the DOS. The above reasoning is 
independent of the probability distribution of A(a;). However, it should be noted that 
A(a;) must not be dominated by one sign. If (p{x) « nvr and (~l)"A(a;) is negative, 
f{x) will fluctuate around the stable position near mr + uj/Ag. ^p{x) = (n + l)7r 
can only be reached if (— l)"A(x) is positive on average over a finite interval. We 
therefore conclude that we expect a Dyson singularity if A(a;) is real and fluctuates 
around Aav = (A(x)) with Aav sufficiently small. 

For complex A (a;), fluctuations of the phase of A (a;) can be mapped via the gauge 



transformation (130) onto a forward scattering potential. Since the amplitude |A(a;)| 
is always positive and the phase fluctuations lead to an effective local shift of the 
frequency w, there should be no Dyson singularity. Instead, we expect a suppression 
of the DOS, i.e. a pseudogap. 

5.2 Numerical algorithm 

In the following, we present an exact algorithm which for stepwise constant potentials 
allows for simultaneous numerical calculations of the integrated DOS and the inverse 
localization length. By choosing the step size sufficiently small, the integrated DOS 
and the inverse localization length may be calculated for arbitrary given potentials. Let 
us partition the interval (0, L) into N intervals (x„,x„+i) of length (5„ — Xn+i — x„ 
with xq ^ < xi < . . . < xn = L, such that A„ = A(a;) and Vn = V{x) for 
x„ < X < Xn+i- Let us also define a)„ as a)„ = w — 

To find an exact anal ytic soluti on f or the given stepwise constant potentials of 



the equations of motion ( |121| ) and ( |122D , let us consider again the related S'-matrix 
iS'(L,0;w), which can be written as the finite product 

N-l N-l 

S{L, 0;u;)=Y[ 5„ = [| 5(a;„+i, x„; w) , (212) 

n=0 n=0 

where 5„ is given for < |a„|2 by Eq. (|l|) and for luI > |A„|2 by Eq. (|2|). Eq. 
( pl2[ ) implies the recurrence relation S{x„+i,0; uj) = 5„ S{x„, 0; uj), which can be cast 
into the following recurrence relations for ipn = V?(x„) and C„ = (i^n)'- 

V9„+i = - 2Im(lnz„) , (213) 
Cn+i = C« + 2Re(lnz„) , (214) 



L. Bartosch, Fluctuation effects in disordered Peierls systems 



45 



where z„ = (571)22 + ('S'n)2i exp(i(^„). Note that Eqs. (213) and (214) are integrated 
forms of the equations of motion ( |121[ ) and (122). The real and imaginary part of lnz„ 
can be determined by writing In z„ as In z„ = In |z„| + i arg (z„). The argument of z„, 
arg (2,1) can be obtained up to a multiple of 27r from 

/ Re z„ ' 



arg (z„) — lunin + sgn [Im2:„] arccos 



To find the integer 
arg (z„ 



27r 



int 



(215) 



(216) 



z„(x) is an analytic function of 
2, it follows from Eq. (|l|) that 



we define Zn(x) by z„ with (5„ replaced by a; — a:„ 
X and ai X — Xn+i agrees with z„. For Q'^ < |A 

Im [z„(a;)] oc sinh[-\/|A„p — oj^{x — x^)] does not change its sign for any x > a;„, such 
that I arg [zn(a;)] | < ti" and m„ has to be zero. For cj^ > | A„p, however, Im [z„(a;)] oc 



5m[^/^}l^^\Kn^ (x-Xn)], such that |[arg(z„(a;)) /7r]j^t| 



|A„|2 {x - a;„)/7r 



Since the constant of proportionality is negative for uj„ > |A„| and positive for 



mt 



|A„|, it follows 



sgn((i„)v^ 



r,2 



27r 



(217) 



int 



To summarize, we can simultaneously calculate the integrated DOS and the inverse 
localization length for arbitrary stepwise constant potentials using the following iter- 
ative algorithm with the initial values ipo = Co = 0, 







- 2 


, FT 1 /Rez„Y 
Zirnin + sgn [Im z„\ arccos -j — j- 1 




Cn+1 


= Cn^ 


f 21n|z„| , 



where z„ and m„ are given by 



Zn = (S'„)22 + (5'«)21 exp (iipn) , 

1 sgn((I'„)^J)2 - |A„|2^^ 



(218) 
(219) 

(220) 
(221) 



and the matrix elements of Sn are determined by Eqs. (|91|) and (|92[). 
5.2.1 Generation of disorder 

In the case of finite correlation lengths, specific extensive physical quantities (which 
are obtained by relating extensive quantities to the length of the system) show a self- 
averaging effect as the length of the chain increases |20|, i.e. they become independent 



27r 



- mt 







|A„|2 >0 
|A„P<0 
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of the concrete reahzation of the disorder. This self-averaging effect can be understood 
by partitioning a very long macroscopic chain into a large number of chains, each one 
still being much longer than the correlation length and any other microscopic length 
scale. In this case, boundary effects between adjacent parts of the original chain 
may be neglected and we are practically left with an ensemble of a large number of 
independent chains. Physical quantities can now be calculated for each chain indi- 
vidually, assuming independently all possible values with their respective statistical 
weight. Specific extensive quantities are now given by the ensemble-average, giving a 
non-random value in the thermodynamic limit. The DOS and the inverse localization 
length can therefore be calculated by generating one typical very long chain. 

Gaussian disorder 

To generate Gaussian disorder at the sample points Xn with the first two moments 
satisfying 

(A(x))=Aav , (A(x)A(x')) = A2e-l^-^'l/« , (222) 

where A(a;) = A.{x) — Aav, we use a realization of an Ornstein-Uhlenbeck process 
described in more general form in Ref. Using the Box-MuUer algorithm |7|], we 
generate independent Gaussian random numbers (/„ with (g„) — and (g^) — 1. For 
real A(a;), we generate A„ = A„ — Aav recursively by 



Ao = A^go , A„+i = an An + ^/l -al Asg„+i , (223) 

where a„ = e"'''"'/^. This Markov process leads to a Gaussian random process with 
the desired correlation functions. The Markov property of the algorithm allows us to 



generate the disorder simultaneously with the iteration of the recurrence relations ( 21S ) 
and (|219| ), so that the algorithm presented above practically needs no memory space 
and, in principle, arbitrary long chains can be considered. If we choose |(Sn|/^ <SC 1 (in 
practical calculations we choose |(5n|/^ ~ 0.0001 to 0.05 depending on A^^ and make 
sure that lessening of |(5„|/^ does not change the results), the (integrated) DOS and 
the inverse localization length may be calculated with arbitrary accuracy numerically. 

Of course, the above algorithm can also be used to generate Vn and in the complex 
case, ReA„ and ImA„ can be generated by replacing A^ by As/V2 (see Ref. |2^). 

5.2.2 Results 

First numerical calculations of the DOS of the FGM in the regime of finite correlation 
lengths were done by myself in collaboration with Peter Kopietz using an algorithm 
similar to the one presented here |2^. Simultaneously, Millis and Monien presented 
their data obtained by an exact diagonalization of a lattice regularization of the FGM 
However, these authors did not make any attempts to relate their results to the 
continuous FGM which would have allowed for a more direct comparison with the 
solutions given by Sadovskii 0. 

In contrast to the algorithm described in Ref. |2^, the algorithm presented here 
does not only allow for a numerical calculation of the (integrated) DOS, it is also 
capable of a simultaneous evaluation of the localization length which for finite ^ has 
never been published before. 
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Commensurate case 

In Fig. 1^, we show our numerical results for the DOS pitLi) and inverse localization 
length £~^{u!) for real A(a;) (with Aav — and V{x) — 0), which refers to the 
symmetric phase of a commensurate system with no forward scattering. Except for 
As^ = 1000, 0.2 we have chosen the same values of the dimensionless parameter As^ 
as m Fig. 7 of Ref. 0. One clearly sees the Dyson singularity in the DOS which exists 
for any finite value of ^ and overshadows the pseudogap at sufficiently small energies. 
One can also see that this Dyson singularity is accompanied by a singularity in the 
inverse localization length. The inverse localization length drops to zero at w = 0, in 
accordance with the exact result £^^{0) — (Aav) [see Eq. ( |144| )]. 

The Dyson singularity in the DOS is missed by Sadovskii's algorithm ||l^. For a 
more quantitative description of the Dyson singularity we have plotted the logarithm 
of the integrated DOS M/Ag versus the logarithm of — ln(ci;/As). For frequencies 
between lu = 10~^^As and ui = lO^^A^, we find that the data can be very well fitted 
by a straight line, such that 

■^'"'-"° ii„(t;S'« ' <^^'" 

which implies for the DOS 

^^"^ = (./A.)|lnti/A.)|i+"(0 ' (225) 
with A{^) — a{^) -B(C)- Pfots of the exponent a(^) and the weight factors of the Dyson 



singularity A(^) and -B(^) are shown in Figs. |10| and 11. For AsS, <C 1 our data is 
consistent with the white noise result a — 2. As AgS, increases, the exponent a{£,) 
decreases, assuming for large correlation lengths ^ the finite value 

a(0 ~ 0.41 , A^ > 500 . (226) 

Fitting the data for A{^) in the regime between AgS, — 500 and Ag^ — 10000 to a 
power law shows that the weight of the singularity of the DOS vanishes as 

yl(e) = 0.175 (A4)^°-^5 ^ (227) 

The plot of the DOS given in Fig. || shows that for large correlation lengths ^ the Dyson 
singularity only overshadows a pseudogap, such that p{uj) takes a minimal value at a 
certain frequency (jJ*{£,)- A double-logarithmic plot of p{uj*) versus (Ag^)"-^ is given 
by the triangles in Fig. |l^. The straight line gives a fit to a power-law: 

p{u:*)/po = C{AgO-^ (228) 
We find 

C = 0.482 ±0.010 , = 0.3526 ± 0.0043 (229) 

The circles in Fig. |l^ show oj* where p{uj) is minimal. The long solid line is a fit to a 
power-law 



u;*/Ag^D{AsO-'' 



(230) 
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Fig. 9 Plot of the DOS p{io) and the inverse localization length £~^{uj) for real A(x) with 
Gaussian statistics, A^L — 10*, and finite correlation lengths Ag^ = 1000, 100, 10, 2.0, 1.0, 0.5, 
and 0.2. For any finite ^, we find p(0) — oo and £~^{lj) = 0. The dashed line represents 
the exact result derived in Section H and for lo > 0.2As is almost indistinguishable from the 
result for A4 = 1000. 



L. Bartosch, Fluctuation effects in disordered Peierls systems 



49 




0.001 0.01 0.10 1.00 10.00 

Fig. 10 Plo t of the exponent a(^) defined 
by Eq. (224) for frequencies between ui = 
10~^As and u = 10""As. While our re- 
sults for small As^ are consistent with the 
white noise result a = 2, in the opposite limit 
A4 > 1 we find Q f» 0.41. 




0.001 0.01 0.10 1.00 10.00 

Fig. 11 Plot of the weight factors A{£) 
and -8(0 of the Dys on singularity defined by 
Eqs. (224) and (225) for frequencies between 
uj = 10^ As and uj = 10""As. As the cor- 
relation length increases, both weight factors 
scale to zero. 



Here, we find 

D = 0.2931 ± 0.0074 , 7 = 0.3513 ± 0.0051 . (231) 

such that within numerical accuracy /i = 7. The proportionality of p{ij->*) to the energy 
scale w*, which can be interpreted as the width of the Dyson singularity, can also 
directly be seen in Fig. |^. Finally we note that for A^^ < 0.2 our algorithm produces 
results consistent with the white noise limit A^^ ^ 1. From the exact solution of 
Ovchinnikov and Erikhman ||l|| we obtain ^(w^J/po 0.9636 and uj* 1.2514 A^^ 
which determines the short solid line in Fig. |l^, describing uj* (^) in the white-noise 
limit. 

Incommensurate case 

The DOS p{uj) and inverse localization length £^^{uj) for complex A(a::) (and Aav = 
V{x) — 0), which refers to the symmetric phase of the commensurate case with no for- 
ward scattering are presented in Fig. Neither the DOS nor the inverse localization 
length involve a singularity. In fact, a direct comparison of our results for the DOS 
with those obtained from Sadovskii's algorithm shows a good agreement. 




1/A4 



Fig. 12 Double-logarithmic plot of 
p{uj*)/po as a function of l/As^ for 
real A{x) (triangles) and complex 
A(a;) (diamonds), where ui* is the en- 
ergy for which the DOS assumes its 
minimum. While uj* = for com- 
plex A(a;), the circles give the double- 
logarithmic plot of u}*/As for real 
A(x) as a function of 
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Fig. 13 Plot of the DOS p{io) and the inverse locahzation length £^^{lu) for com- 
plex A(a;) with Gaussian statistics, AgL = 10*, and finite correlation lengths A^^ — 
1000,100,10,2.0,1.0,0.5, and 0.2. The dashed line represents the exact result derived in 
Section ^. This line is hardly recognizable because it is almost indistinguishable from the line 
for A4 = 1000. 
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For a more quantitative comparison, the diamonds in Fig. show the DOS p(0) 
at the Fermi energy. A fit to a power-law gives 

p(0)/po = C(A,0"^ (232) 

with 

C ^ 0.6397 ± 0.0066 , /i = 0.6397 ± 0.0024 . (233) 

Note that within numerical accuracy we find C = ^. This result should be compared 
with Sadovskii's approximate result C = 0.541 ± 0.013 and /i = 1/2. 



5.3 Phase fluctuations only 

As already discussed at the end of Section |[ far below the mean-field critical tempera- 
ture , amplitude fluctuations of the complex order parameter A(a::) = |A(x)|e"'(^) 
are frozen out and only phase fluctuations survive. In terms of the "superfluid velocity" 
V{x) — dxi}{x)/2, the process of averaging can be written as 



e 



-PF{V} 



JV{V}e'Pnv} 



(234) 



where is the free energy functional given in Eq. 1^3^). Since -^{1^} is Gaussian 

and local, the process of averaging is described by Gaussian white noise. The first two 
moments of V{x) are given by {V{x)) = and 

{V{x)V{x'))^2^^5{x-x'), (235) 

with [see Eq.|9) £,{T) = sps{T)/2T. WeU below the mean-field temperature T^'^ , we 
can use the BCS gap equation As ~ 1.764 T,?^^ and Ps{T) w po = ""^^ to get for the 
dimensionless parameter^ As£,{T) 

AsC{T) = 0.281 s T^^/T . (236) 



As already shown in Section 0, Eq. (235) implies the exponentially decaying correlation 



function {A{x)A*{x')) = Af exp (-|j; - x'|/$(T)). 

To calculate physical quantities like the DOS or the inverse localization length, we 



use the gauge invariance of these quantities under the gauge transformation (130) and 
map the phase fluctuations of the order parameter A{x) = Age"'*^^-* onto the effective 
forward scattering potential V{x) — dx'd{x)/2. 

5.3.1 Density of states and inverse localization lengtli 

The DOS and the inverse localization length for the remaining problem involving only 
a constant gap parameter and forward scattering described by Gaussian white noise 
were already calculated in the previous section. With the exception of the constant 
shift in the inverse localization length, the results are identical with those for the 
incommensurate case without forward scattering and (A(x)) = Aq ^ 0. Substituting 

-"^^The prefactor can be expressed in terms of the Euler constant 7, such that As^(T) = sT^^ /2e'^T. 
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in Eqs. (|T9|) and ( pOOj ) Aq by A^, £> by 1/4^(T) and setting D = (since there is only 
forward scattering) which imphes a completely different interpretation of the resulting 
equations, we get T{uj) — —iuj + As Ii-i4Lui (4As^)//_j4^^ (4As^), or, equivalently, 



It follows from Eq. (201) that the integrated DOS is given by 



sinh {Attu!(,) 
4< 



|/,4^^ (4A,0r 



(237) 



(238) 



Plots of the DOS and the inverse localization length for characteristic values of A^^ 
are shown in Fig. ^ At the Fermi energy, the DOS simplifies to 
Po 



PiO) 



(239) 



such that, as the temperature is lowered and the correlation length grows, the DOS 
at the Fermi energy vanishes exponentially, 



p(0) - SnpoA^ exp(-8A,0 , 4A„e > 1 



(240) 



This result is in contrast to the power-law behavior of the DOS as predicted by Gaus- 
sian statistics. A plot of the DOS at the Fermi energy is shown in Fig. [ij as a function 
of 1/AsS,. For a comparison, we have also plotted p(0) for Gaussian statistics and the 
result found in the Born approximation, which at low temperatures can only poorly de- 
scribe the quantitative behavior of the pseudogap. For T — T^^/4 which corresponds 
to As£, ~ 2.0, we find that the DOS p{Q) for phase fluctuations only is less than 10~'^po 
while the Born approximation and the numerical exact result for Gaussian statistics 
suggest a value of order po/4. 

As can also be seen in Fig. there is no pseudogap for A^^ ^0.1. Note, however, 
that these correlation lengths correspond to temperatures of order , where ampli- 
tude fluctuations are important. Nevertheless, using only phase fluctuations gives a 
good qualitative description of the DOS for all temperatures T. 

For temperatures well below the mean-field temp erat ure where our theory 

becomes quite accurate, we have, according to Eq. (236), A^^ ^ 1, such that the 
Bessel function liyivz) — e^'^'^/^ Jii,{ivz) may be approximated by an Airy function. 
For large correlation lengths we expand the resulting equation for oj around As and 
obtain to leading order in 1/4 A^^ a maximum of the DOS piyi) at As, described by 
the inverted parabola 



p{uj) - Po 



a(4AsO'/'-M4A.O 



5/3 



a: 



with a 



1 C2 



0.731 and b 



1 1 



1 



- 1 



(241) 



0.258, where ci = Ai (0) = 



3-2/3/p(2/3) ~ 0.355 and C2 = - Ai' (0) = 3-i/Vr(l/3) w 0.259. Note that Eq. ( |24l] ) 
implies that the maxinram of the DOS diverges as (4As^)^/'^ oc T~^^^. 
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Away from w = Ag, we find for u> > Ag 



/\2 

p{u) ~ Po ( 1 + ^ 



3/2 

1 I > 1 , 



(242) 
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Fig. 15 Plot of the DOS p(0)/po 
as a function of 1/As^. The solid 
line gives the DOS for phase fluctu- 
ations only while the dashed line is 
the result found in the leading order 
Born approximation [see Eq. pGl and 
the diamonds give the DOS evalu- 
ated for Gaussian statistics [see Sub- 



section 5.2 [ 



which is independent of A^^ and agrees with the mean-field result. And for w < 
we find 



p{uj) ^ %pq\J /S.1 — ^ arccos (w/ As) [1 -I- exp(— Stto;^)] 

X exp ^— 8[\/A2 — a;2 ^ — arccos (tj/A^)]^ , 

/ X 3/2 

4A,C (l - ^ j » 1 • (243) 

If w <C As and 4a;^^/As ^ 1, this result simplifies to 

p(0) - Sttpo As^ cosh(47rtjO exp(-8AsO . (244) 
5.3.2 Pauli paramagnetic susceptibility 

The Pauli paramagnetic susceptibility is defined as the contribution of the conduction 
electrons to the susceptibility and can be written in terms of the DOS. A magnetic 
field H shifts the energy levels of the electrons by an amount ±fiBH, where /is is 
the Bohr magneton and the sign depends on the spin orientation of the electron with 
respect to the field. The resulting different occupation of spin-up and spin-down states 
leads to a magnetization density M{T) which for small magnetic fields is linear in H . 
If p{ijj) = p{-ijj) as in the FGM, the susceptibility x{T) = dM{T)/dH can be written 
as 

Placing the asymptotic expression ( 243 ) with ^(T) = Ps/T into this equation, we find 
r , n3 ? ^ 1 + exp[-8^psw/r] 



x{T)Ixq - ^ duj VA2 - tj2 arccos 



A 

X exp 



T 



-— 1 — 8ps arccos — 
As V V As 



(l-fexp[-2cj/r])2 
+8ps^l-{u;/Asf ) , (246) 



where xo = 2^^jOo- For ps > po/4 (and T <^ T^^), the integrand is sharply peaked 
at oj = cos(l/8ps)As, such that the integral may be evaluated by a saddle point 
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Fig. 16 Plot of the susceptibility 
X(r) calculated for ^{T) = Ps{T)/T 
[with Ps{T) given by Eq. (p7j)] and 
As(T) determined by the minimum 
of the Ginzburg-Landau functional. 
For a comparison, we also show as the 
dashed line the DOS at the Fermi en- 
ergy, p(0). 



integration, resulting in 



X{T) nr /sin(l/8p,)A,\^/' / 8p3sin(l/8p3)A, 



exp 



Xo \ Ps \ T J \ T 

T « 4p,A, (1 - cos(l/8p,))'/' 



(247) 



Note, that the temperature restriction is necessary for the asymptotic expansion of 
the DOS to be valid. Although the susceptibility x{T) vanishes exponentially, the 
exponent 8ps sin(l/8ps)As/r is smaller than the exponent SpsA^/T, which governs 
the DOS at the Fermi energy. However, as ps approaches po/4, the two exponents 
become identical, such that for ps < po/4 the DOS and the susceptibility have the 
same exponential dependence on T. 

A numerical evaluation of the susceptibility xiT)/xo for ^(T) — ps(T)/T given by 
Eq. ( p7| ) and As(T) determined by the BCS gap equation ( ^8| ) is shown in Fig. ^ 
For a comparison, we have also plotted the DOS p{0)/po as a function of temperature. 
The two are not identical because for small temperatures (and ps = po), the major 



contribution to the integral in Eq. (245) comes from the frequency region just below Ag. 
In Fig. we also show susceptibility data taken from Ref. [ p2| for incommensurate 
quasi one-dimensional conductors which undergo a Pererls transition. 

It is quite surprising that our plot of the susceptibility is very similar to the plot 
obtained by Lee, Rice and Anderson ^ which perfectly fits experimental data j|, 
^,|8^. However, to explain experimental data, Lee, Rice and Anderson ^ had to 
base their calculations on a real order parameter with a correlation length which 
for low temperatures increases exponentially as the temperature is lowered. Only 
the exponentially increasing correlation length of a real order parameter could lead 
to an exponentially decreasing susceptibility and the prediction of Tj^^ « T^^/4. 
Here, we have shown that these predictions should also hold for a complex order 
parameter with a correlation length which increases as 1/T. Since most Peierls chains 
are incommensurate and the susceptibility of many incommensurate Peierls chains 
has been compared with the theory by Lee, Rice and Anderson our results are 
of major experimental relevance. For a comparison between theory and experiment, 
it should be recalled that we have only used a strictly one-dimensional model with 
phase fluctuations only. At higher temperatures, one should also include amplitude 
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fluctuations. 



5.3.3 Thermodynamic quantities 

The DOS encapsulates the whole thermodynamics. Let us first consider the electronic 
free energy Fei{T) with respect to the gapped state for which we have Poo(^) = 
po9{Lo^-Al)\Lo\/{Lo^-Aly. 



Partial integration leads to 



sL 

T 



(248) 



1 



g/3c^ + 1 



spoLIm J did [T{uj) - T^{uj)] ^ ^ , 



(249) 



where r(a;) = £~^{u}) — inj\f{uj) and roo(ij-') = y/ A'^ — {cu + iO)^. The square root has 
to be taken such that Too{oj) —iuj for co oo. Since r(a;) is analytic in the upper 
half plane, the integral may be done by closing the integral in the upper half plane 
and using the residue theorem. We find 



F,,iT)~Ft°^{T) = -spoL^ J2 [Rer(zi„) 



(250) 



For the FGM with phase fluctuations only, we obtain by placing Eq. (237) into this 
equation 



F,i{T)-Ft^{T) = spoLA 



27r ^ 

w„>0 



r' 



Up to an irrelevant additive constant, F^^ °°(T) is given by 



(T) - -spoL2Al 



du yi. 



1 



1 



,1 eA,n/T^.l' 

which for small temperatures is exponentially small: 



(251) 



(252) 



(253) 



While in the general case we have to add Eqs. (251) and (252) to get the free energy 
i<"oi(T), for low temperatures we can neglect the exponentially small contribution given 
by Eq. ( ^53[ ) , such that F^\ (T) is determined by the right-hand side of Eq. ( |25lD . 

For T <C T^^, the dimensionless correlation length A^^ is large, and a uniform 
asymptotic expansion of /^(j/z) and Il{yz) [see A&S, Eqs. (9.7.7) and (9.7.9)] can be 
used to flnd for the leading terms of the free energy 



F,,{T) ^ spoLAl 



1 



1 



1 



4 4A,e 12 (4A,02J 



(254) 
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Electronic speciGc heat 

An experimentally accessible thermodynamic quantity is the electronic specific heat 
which can be expressed in terms of the free energy as 



Cei(r) = -T 



(255) 



The low-temperature behavior of Coi(T) can be obtained from Eq. (254): Using ^(T) 
s/9s(0)/2r, it directly follows 




(256) 



where the specific heat of free electrons is given by 



C^AT) = s — p^LT 



(257) 



Although the DOS exhibits a pseudogap and vanishes exponentially near the Fermi 
energy as the temperature is lowered, the electronic specific heat Cc\{T) vanishes only 
linearly in T, as for free electrons. 

We conclude this section with a summary of the central results of the FGM valid 
at low temperatures where phase fluctuations dominate in Table |^. 



Table 1 Asymptotic low temperature results for the FGM describing electrons with spin. 
Note that we have reintroduced the Fermi velocity vf and note also that the mean-field 
critical temperature serves as the only energy scale. For generalizations of the formulas 
see the text. 



superfluid density ps (T) ^ po 



ttvf 



correlation length 



VF 

ttT 



, . , p(0) U.IT^^ ( 4.49 T,MF 
density oi states ~ — exp ' 



Po T '\ T 

1.76 T^^F 



inverse localization length £ ^(0) 

susceptibihty - 1.74 ( ) exp 

electronic specific heat 



X{T) ^JT^^V'^ ( l.l2Tf^ 



Cei(r) 1 



CO(r) 32 
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6 Conclusion 

In this work, we have discussed the density of states (DOS) of the fluctuating gap 
model (FGM) and related quantities like the inverse localization length, the Pauli 
paramagnetic susceptibility and the low-temperature specific heat. We introduced the 
FGM as an effective low-energy model describing the electronic properties of Peierls 
chains and emphasized the fact that the FGM also finds its applications in other 
physical contexts: Spin chains can be mapped by a Jordan- Wigner transformation 
onto the FGM and in order to explain the pseudogap-phenomenon in underdoped 
cuprates above a phase transition, higher-dimensional generalizations of the FGM 
have been used. 

With the rediscovery of the FGM in the context of high-temperature superconduc- 
tivity, a previously unnoticed subtle error in Sadovskii's widely used Green function 
of the FGM was brought to light. This error called for a reinvestigation of the FGM. 

After setting up a non-perturbative theory which, in principle, allows to express 
the one-particle Green function as a functional of an arbitrary given realization of the 
disorder, we derived a simple equation of motion whose solution determines the DOS 
and the inverse localization length. Starting from this equation, we could rederive all 
known results for the FGM in the white noise limit. 

Considering the equation of motion governed by the phase which determines the 
DOS, we argued that the Dyson singularity found in the white noise limit for com- 
mensurate Peierls chains should not be an artifact of the white noise limit, but should 
be present for any finite correlation length in contradiction to Sadovskii's solution. 
Our following numerical calculation of the DOS and inverse localization length con- 
firmed this prediction and showed also that for large correlation lengths, the Dyson 
singularity only overshadows a pseudogap. Although Sadovskii's algorithm misses this 
singularity, his solutions for the incommensurate case where there are no singularities 
in the DOS give a fairly good approximation to the exact result. 

In the pseudogap-regime below the mean-field critical temperature, fluctuations of 
the order parameter cannot be described by Gaussian statistics. Instead, as the tem- 
perature is lowered, amplitude fluctuations get gradually frozen out, and the amplitude 
takes on a value given by the minimum of the Ginzburg-Landau functional and only 
long-wavelength gapless phase fluctuations survive. Using a gauge transformation to 
map the phase fluctuations of the order parameter onto an effective forward scatter- 
ing potential, we could even find an exact solution for the FGM involving only phase 
fluctuations which should be valid in the low temperature regime. We found that the 
low-tcmperature specific heat is linear in T and that both the DOS at the Fermi energy 
and the Pauli paramagnetic susceptibility vanish exponentially as the temperature T 
is lowered, the ratio of the former to the latter also vanishing exponentially. The Pauli 
paramagnetic susceptibility has been measured in various experiments and is in good 
agreement with our results. 

Having discussed quantities related to the DOS, one would also like to calculate 
quantities like the spectral function. This has been done for a special non-Gaussian 
probability distribution involving amplitude and phase fluctuations in Ref. [ p2[ , but 
accurate results for realistic probability distributions (e.g. for phase fluctuations only) 
are not known yet. 
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